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0)  =  scalar  ratio  following  from  the  vanishing  of  the  Jacobian  determinant. 

=  infimum,  or  lower  bound  of  characteristics  in  a  simple  wave. 

6’  =  control  volume  enclosing 

T  =  mapping  from  x-t  space  to  T^~  U  space,  where  T  is  the  local  stagnation 
temperature. 

{L  ]  =  set  of  lines  in  domain  3)  of  function  /  . 

a  a=l,  2,  ... 

0  =  origin  of  coordinates  in  x-t  space. 

c'(y)  =  derivative  of  differentiable  function  C(y)  with  respect  to  y. 

AOj  =  interval  between  two  points  on  a  characteristic  of  type  I. 

K  =  ratio  of  the  distance  of  a  given  point  in  interval  Ao^  from  the  left-most  point  in  the 
interval  to  the  length  of  the  interval. 

0(Ao^)"  =  n’  th  order  error  on  interval  Aa^  . 

SB  =  control  volume  having  the  same  cross-section  as  the  void  main,  extended  to  the  left, 
into  the  infinite  reservoir. 

L  =  length  of  control  volume  SB  extending  to  the  left,  into  the  infinite  reservoir  (See 
Appendix  for  this  usage). 

L  ^  =  length  of  control  volume  SB  extending  to  -<»  . 

o  =  right-most  face  of  control  volume  SB  ,  which  is  the  surface  in  which  the  flow  is 
choked. 

=  remainder  of  surface  of  control  volume  SB  when  is  excluded. 

a  =  entire  surface  of  control  volume  SB  (See  Appendix  for  this  usage). 
r  -  ratio  of  absolute  pressure  in  the  quiescent  flow  region  in  the  reservoir  to  the  absolute 
static  pressure  at  the  plane  of  choked  flow  (See  Appendix  for  usage). 

B  =  ratio  of  absolute  temperature  in  the  quiescent  flow  region  in  the  reservoir  to  the 

absolute  static  temperature  at  the  plane  of  choked  flow  (See  Appendix  for  usage). 
P(z)  =  function  of  pressure  ratio  (See  Appendix  for  usage). 
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1. 


SUMMARY 


Temperature  data  acquired  from  turbine  stator  transient  flow  tests  performed  in  the 
Turbine  Research  Facility,  operated  by  the  Propulsion  Directorate  of  the  Air  Force  Research 
Laboratory,  show  a  large  initial  thermal  impulse  appearing  in  each  test.  This  impulse  is 
clearly  not  related  to  the  performance  characteristics  of  a  particular  test  article,  and  it  raises  a 
question  concerning  its  interpretation  with  regard  to  the  test  data  in  which  it  is  imbedded. 

The  flow  physics  producing  the  thermal  impulse  is  assumed  to  be  associated  with  the 
sudden  expansion  of  the  working  fluid  from  a  large  reservoir  into  the  domain  of  the  test 
article,  which  is  initially  in  vacuo.  In  August  of  1997,  a  one-dimensional  mathematical 
analysis  of  the  wave  phenomena  associated  with  the  expansion  of  a  perfect  gas  from  an 
infinite  reservoir  into  a  void  main  was  undertaken  in  order  to  understand  the  origin  of  this 
impulse,  or  spike;  and  to  understand  the  variables  upon  which  it  depends.  It  was  understood 
at  the  outset  that  this  model  would  apply  only  to  the  very  early  development  of  the  flow  in 
the  actual  main  supplying  the  test  article  with  air  from  the  facility’s  reservoir. 

The  gas  flow  from  the  infinite  reservoir  into  the  main  was  found  to  be  embodied  by  a 
simple,  centered,  backward-facing  expansion  wave.  The  limit  of  characteristics  bounding  the 
wave  at  its  tail  was  foimd  to  be  a  shock,  since  the  Rankine-Hugoniot  Relations  are  satisfied  in 
a  non-trivial  way  across  that  bound.  The  lead  characteristic  of  die  expansion  wave  was  found 
to  be  permanently  located  at  the  exit  plane  of  the  reservoir  following  initiation  of  the  flow 
process.  In  the  one-dimensional  flow  model,  it  is  kept  in  place  by  an  isentropic  expansion 
shock  which,  in  the  case  of  an  infinite  reservoir,  does  not  contradict  the  Second  Law  of 
Thermodynamics,  and  thus  is  physically  permissible.  This  assertion  is  demonstrated  in  the 
appendix  of  this  report.  The  stationary  position  of  the  lead  characteristic  of  the  backward¬ 
facing  expansion  wave  implies  that  flow  velocity  is  precisely  critical  at  initial  time,  and  for 
any  time  interval  following  initial  time,  since  the  main  is  infinite  in  the  mathematical  model. 

The  stagnation  temperature  was  found  to  be  a  strict  function  of  soimd  speed 
throughout  the  domain  of  the  wave,  and  was  found  to  increase  monotonically  from  a 
minimum  at  the  critical  point,  at  the  head  of  the  wave,  to  a  maximum  at  the  tail  of  the  wave. 
The  ratio  of  the  maximum  to  the  minimum  value  was  found  to  be  fixed,  and  to  have  a  value 
of  six  for  standard  air.  The  strict  dependence  of  stagnation  temperature  upon  sound  speed  was 
found  to  be  due  to  the  fact  that  the  Riemann  Invariant  on  each  spanning  characteristic  of  the 
simple  wave  has  precisely  the  same  value.  The  stagnation  temperature  is  fixed,  but  has  a 
distinct  value  on  each  straight  characteristic  of  the  simple  wave. 

It  is  concluded,  based  on  the  findings  of  this  analysis,  that  the  observed  thermal  spike 
occurring  in  the  blow-down  process  is  the  result  of  the  heat  of  the  compression  to  form  the 
shock  proceeding  into  the  main  when  the  flow  is  initiated.  The  calculated  stagnation 
temperature  distribution  in  the  rarefaction  wave  is  the  artifact  of  that  process.  The  theory  of 
characteristics  does  not  yield  any  information  on  the  flow  physics  associated  with  the  shock, 
itself,  but  it  accurately  portrays  the  compression  process  everywhere  in  the  domain  of  the 
flow  following  the  shock. 
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2. 


INTRODUCTION 


The  Turbine  Research  Facility,  operated  by  the  Propulsion  Directorate  of  the  Air  Force 
Research  Laboratory,  is  a  transient  blow-down  facility  which  has  the  purpose  of  obtaining 
accurate  aerodynamic,  thermodynamic,  and  thermal  performance  data  from  large  turbine 
stages  without  the  expenditure  of  impractical  amoimts  of  power.  This  objective  is  clearly  met 
by  this  facility,  and  it  is  not  the  intention  of  this  report  to  discuss  its  configuration  or 
operation.  This  is  well  documented  in  open  literature.  However,  an  interesting  thermodynamic 
phenomenon,  sometimes  called  the  "thermal  spike",  appears  at  the  initiation  of  the  blow¬ 
down  process.  The  purpose  of  this  report  is  to  discuss  the  mathematical  physics  of  the 
thermal  spike  in  order  to  better  understand  its  origin. 

The  general  hypothesis  presented  in  this  report  is  that  the  thermal  spike  is  produced  by 
a  backward-facing  rarefaction  wave  which  has,  at  its  tail,  the  initial  shock  proceeding  into  the 
domain  of  the  test  article  when  the  blow-down  process  is  initiated.  The  mathematical  model 
presented  here  is  that  of  a  one  dimensional,  centered,  backward-facing  rarefaction  wave 
created  by  the  expansion  of  a  calorically  perfect  gas  from  a  reservoir  of  infinite  extent  into  a 
void  main.  The  gas  flow  is  assumed  to  be  inviscid. 


3.  THE  FIRST  AND  SECOND  LAWS  OF  THERMODYNAMICS  AND  THE 
DERIVATION  OF  THE  EQUATIONS  OF  FLUID  MOTION 


Consider  Figure  1. 
Here  is  shown  an  infinite 
reservoir  from  which  a  shock, 
followed  by  a  rarefaction, 
issues  into  a  main  in  vacuo  at 
time  0.  The  flow  is  assumed 
to  be  one-dimensional,  and 
perpetually  choked  at  the 
indicated  plane  0-0  in  Figure 
1  after  time  0.  Since  the 
flow  is  one-dimensional, 
there  can  be  no  convergence 
of  the  streamlines  as  the  flow 
enters  the  rarefaction;  so  the 
flow  cannot  accelerate  on  a 
gradient  as  it  approaches  the 
plane  of  choked  flow.  It  is 
therefore  quiescent  to  the  left 
of  plane  0-0,  mamtaining  a 
time-steady  stagnation 


Figure  1 
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thennodynamic  state,  defined  by  fi[xed  stagnation  temperature  and  pressure.  It  turns  out  that 
essentially  no  inconsistency^  is  produced  by  assuming  an  abrupt  acceleration  to  critical  flow 
at  plane  0-0  in  Figure  1,  followed  by  supersonic  flow  throughout  the  rarefaction.  The  thermal 
spike  is  evidently  produced  by  the  flow  work  done  on  the  various  elements  of  the  expansion 
wave,  also  called  a  rarefaction  wave,  which  propagate  at  different  speeds  in  the  rarefied 
region. 

The  analysis  of  the  wave  propagation  problem  requires  that  we  examine  the  First  Law 
of  Thermodynamics  in  its  convective  form.  We  begin  by  considering  control  volmne  V  shown 

in  Figure  2.  Over  any  given  time  interval,  bt  =  ,  beginning  at  state  1  and  ending  at 

state  2,  the  First  Law  of  Thermodynamics  may  be  stated  according  to: 


E 

in 


+  Q  =  W 

CV  CV 


out 


(1) 


where  W  represents  the 

CV  ^ 

work  in  all  forms  done  by  the 
system  upon  its  environment. 

We  may  assume  that  our  system 

is  adiabatic,  and  that  W  is 

strictly  flow  work.  Let 
6m  =  (m,  -  m.)  ,  and 

in  2  1  in 

6m  =  (m^  -  m.)  .  Clearly, 

out  2  I' out 
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- - - 

■  E»  tpil 

in 

E  out 

- > 

1  1  1  CV 
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> 

f 

W  CV  O  CV 


(m  L  -  (m  ),  =  6m.  -  6m„  . 

'  cv'2  '  cv-'l  in  out 


Control  Volume  V 
Figure  2 


In  this  case,  over  time  interval  6f  ,  it  may  be  written  that: 


W  =  -  p  'V  'bm  +  p  ’V  ’bm 

CV  in  in  in  out  out  out 


(2) 


Note  that  the  specific  total  energy  is  defined  as: 

(V-f) 
e  =  u  +  - 

2 


(3) 


'  Refer  to  Appendix. 
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Keeping  Equation  (3)  in  mind,  we  should  write: 


(^2  -  ^l) 


cv 


6(m  -e  ) 

'■  cv  «- 


(4) 


where  m  and  e  are  the  mean  values  in  control  volume  V. 

CV  CV 

Also: 


E  =  6in.  ‘C 

in  in 


E  =  6ot  •  e  , 

out  our  our 


(5.i) 

(5.ii) 


In  view  of  the  adiabatic  hypothesis,  Equation  (1)  implies  that  the  following  rate  equation 
exists: 


5m. 

Ill 

~67 


(e.  +  p.  -v.  )  = 

'‘in  ^  in  m' 


bE 

Cl 

IT 


bm 


out 


bt 


(e  +  p  "V  ) 

^  out  *  oi/f' 


out  out' 


Each  of  the  rates  in  the  above  equation  has  a  limit  as  bt  becomes  infinitesimal,  resulting  in 
the  expression  written  next: 


dE  dm 


out 


dt  dt 


dm. 

in  V  rv 

(e  +  p  -v  )  -  - •  (e.  +  P.  ‘V,  )  =  0 

^  out  ^ out  out'  ^  in  in' 


(6) 


dt 


But  if  a  is  the  surface  area  of  control  volume  V,  it  must  be  that: 


dm  ,  dm. 

out  V  / 

0^  ^^out  ^  P out  ^ou^  0^  ^  Pin  ^ 

pv)da 


(7) 
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The  role  of  the  flow  work  upon  the  rate  of  change  of  total  energy  within  V  is  made  apparent 
by  Equation  (7). 

Consider  point  z  in  control  volume  V  of  Figure  2.  Let  us  suppose  that  every 
neighborhood  a  oi  z  ,  including  control  volume  V,  is  a  regular,  simply-connected  domain 

upon  which  V  is  continuously  differentiable,  and  upon  which'  e  ,  p  ,  and  p  are  of 
class  C'.  Let  us  further  assume  that  z  is  in  an  Eulerian  reference  frame,  with  the  fluid 

elements  flowing  past  at  velocity  V  .Let  de  be  an  infinitesimal  element  of  a,  and  pick 
control  volume  V  itself  as  an  example.  Then  Gauss’  Divergence  Theorem  implies  that  for  V : 

pv)da  =  IJ^(pv-iiyeda  *  lfp-(V-t}do  = 

-  (8) 

fffjy'(pvye  +  (pvyvejde  + 


Equation  (8)  applies  to  any  neighborhood  of  Z  . 

It  is  apparent  that  if  one  picks  another  point  in  V,  such  as  z^  then  every 

neighborhood  a'  of  z  is  also  a  regular,  simply-coimected  domain  upon  which  V  is 

continuously  differentiable,  and  upon  which  e  ,  p,  and  p  are  of  class  C'.  Again,  we  may 
select  control  volume  V  as  the  example,  and  let  us  divide  it  into  a  set  P  of  contiguous,  non¬ 
intersecting  sub- volumes  8^. ,  called  a  partition,  such  that  each  contains  a  point  in  V  in  its 

interior.  It  is  clear  from  the  properties  assigned  to  V  that  for  a  real  number  e  >  0  ,  partition 
P  may  be  constructed  such  that  the  sum  shown  in  Relation  (9),  taken  over  a  refinement  of  P, 
satisfies  the  indicated  condition: 


— V  5(p^) 

y .  - *e 

i=l,a  dt  ‘ 


dt 


<  € 


(9) 


where  a  is  the  number  of  sub-volumes  in  the  refinement.  It  follows  that: 


‘  A  function  on  a  domain  ®  in  Kp  is  of  class  C  if  all  of  its  partial  derivatives  exist  and  are  continuous  on  &. 
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(10) 


dE 

CV 

'W~ 


d(pe) 

dt 


dt 


Introducing  the  substantial  derivative,  Equation  (10)  becomes: 


Cl 

“aT 


fl 


(11) 


But  we  find  that  Equation  (8)  may  also  be  written  in  the  form: 
J'J'  (pV’nyeda  +  JJp‘(V'n)da  = 

f f ^  Iff 


(12) 


Hence,  Equation  (6)  requires  that: 


+  V  (P  V) 


dt 


=  0 


(13) 


Equation  (13)  guarantees  that  at  any  point  z '  chosen  in  V,  the  following  convective^ 
form  of  the  energy  equation  holds: 

p-—  +V'{pV)  =0  (14) 

Dt 


For  suppose,  at  some  point  z '  in  control  volume  V,  there  exists  some  number  k  such  that 
|k|  >0,  and  such  that: 


*  Convective  forms  are  those  containing  the  substantial  derivative. 
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P* -  +7'ipV)  =  K 

Dt 


Then  there  is  a  sub-volume,  e .  c  ]^  ,  and  a  neighborhood  of  type  a  such  that  8^.  c 
and  which  has  the  property  that,  for  some  |  d  |  >0: 

If 

J  J  J  g 

i 

Note  that  imder  the  hypotheses  used  to  formulate  Equation  (13),  control  volume  V  was  picked 
without  loss  of  generality,  i.e.,  it  was  not  assigned  any  special  properties.  It  follows  that,  for 
sub-volume  : 

iff 

J  J  J 

i 


Hence,  our  supposition  leads  to  a  contradiction.  Equation  (14)  characterizes  the  relationship 
between  local  flow  work  done  and  specific  total  energy  at  a  point  in  a  flow  fleld  for  an 
adiabatic  flow  process;  which  is  the  reason  for  its  derivation  in  this  report. 

We  reiterate  that  by  convective  form  of  an  equation,  we  mean  that  it  is  written  in 
terms  of  its  substantial  derivatives.  There  are  four  convective  equations  which  apply  from  the 
physics  principles  governing  the  flow  field  being  considered  here.  They  are  listed  next. 

First  Law  of  Thermodynamics 

£l^lY-ipV)=0  (15.i) 

Dt  p 

Momentum  Equations 

—  +  -Vp  +A  =  0  A  s  0  (ISH) 

Dtp 


p- —  +  V*(pF)  de  =  0 
Dt 


p- —  +  V'ipV)  dt  = 
Dt 
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^  .  pvf  =  0 

Dt 


Mass  Flow  Continuity 


Second  Law  of  Thermodynamics 


Ds 

Du 

Dv 

+  P-  — 

Dt 

Dt 

Dt 

Continuing,  we  note  that  the  above  system  of  equations  can  be  reduced.  Note  that 
Equation  (IS.iiii)  can  be  rewritten  in  the  form: 


Ds  Du  p  Dp 

y*. _  —  _  ^ 

Dt  Dt  p^  Dt 


(16) 


Making  use  of  the  mass  flow  continuity  equation,  this  equation  becomes: 


r 


Ds 

'Wt 


Du 


V 


(17) 


In  view  of  Equation  (3)  the  substantial  derivative  of  the  internal  energy  may  be  written  as: 


Du 

~Di 


De 


(18) 


Taking  account  of  Equation  (15.ii),  the  substantial  derivative  of  the  internal  energy  can  be 
written  in  terms  of  the  momentum  constraint: 


Du 

iDt 


De 

'm 


+  K* 


+  F*A 


(19) 


The  result  of  these  operations  is  that  Equation  (IS.iiii)  becomes: 
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(20) 


r 


Ds 

'Dt 


£i  +  J:  V  ip  V)  +  V'A.  =  7*  A 
Dt  p 


Making  use  of  Gauss’  Divergence  Theorem,  Equation  (20)  implies  that  for  sub-volume, 
e,-  c  V  >  the  entropy  efflux  meets  the  condition  that: 


d(ps)  pV'A 


dt 


de  +  jj  sipV)’nda  =  0 


(21) 


It  is  immediately  apparent  from  the  form  of  Equation  (21)  that  the  only  lost  work  is  due  to 
the  dissipation  force  per  unit  volume,  A  .  It  follows  that  in  case  A  =  0  ,  i.e.,  the  Euler 
Equations  hold,  the  lost  work  function  must  vanish,  and  the  thermodynamic  process  must  be 
reversible.  If  we  suppose  that  this  is  the  case,  then  our  hypothesis  that  the  system  is  adiabatic 
implies  that: 


de 


=  0 


(22) 


Hence,  if  the  Euler  Equations  hold  throughout  control  volume  V,  and  the  adiabatic 
assumption  holdsS  then,  at  every  point,  p  ,  in  V: 


dt 


In  addition.  Equation  (20)^  is  reduced  to: 


Dt 


*  These  relations  apply  in  a  domain  where  pT>0. 
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at  every  point  in  V.  Now  recall  that: 


Ds 

'Wt 


ds 

dt 


+  V'Vs 


Hence,  we  must  conclude  that  the  following  equation  holds  at  every  point  in  V: 

F-Vs  =  0  (23) 


Let  us  add  the  hypothesis  that  the  Euler  Equations  of  Motion  hold  for  our  description 
of  the  rarefaction  shown  in  Figure  1.  It  is  now  clear  that  the  convective  relations  which  apply 
from  the  physics  principles  governing  the  flow  field  are  three  in  number.  These  governing 
equations*  are: 

Momentum  (Euler  Equations) 

+  Ivp  =  0  (24.i) 

Dt  p 


Mass  Flow  Continuity 


—  +  pV'F=0  (24.ii) 

Dt 


Ds 

~Di 


=  0 


ds 

dt 


=  0 


Entropy  Convection 


(24.iii) 


Finally,  we  note  that  for  a  compressible  medium  in  which  pressure  is  strictly 
dependent  upon  density,  sound  speed  is  expressible  as: 


c  = 


dp 


^  dp 


(25) 


‘  These  relations  apply  in  a  domain  where  pT>0. 
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4. 


CHARACTERISTIC  THEORY  APPLIED  TO  THE  ONE-DIMENSIONAL  WAVE 


Letting  the  subscripted  variable  represent  partial  differentiation,  the  inviscid  equations 
of  fluid  motion  described  in  Section  3,  when  reduced  to  one  space  dimension,  appear  as 
follows: 


P,  +  U’P^  +  P’U^  =  0 

(26.i) 

u  +  U'U  +  --P  =  0 

(26.ii) 

f  ^  p  X 

II 

o 

(26.iii) 

O 

II 

1 

JC 

(26.iiii) 

Clearly,  Equations  (26.iii)  and  (26.iiii)  imply  that  the  one-dimensional  expansion  wave 
described  here  is  isentropic,  since  U  >  0  everywhere.  When  Equation  (25)  is  introduced, 
the  above  system  reduces  to  two  equations  in  the  following  form: 

Q  +  U-p  +  p-U  =  0  (27.i) 

X  X 

p-u^  +  P'U'U^  +  c^-p^  =  0  (27.ii) 


We  suppose  that  the 
partial  derivatives,  p^, 

p  ,  U ,  and  U  all  exist, 

X  t 

and  are  continuous 
throughout  domain, 
shown  in  Figure  3.  Then 
function,  / ,  solving  the 
above  equations,  is  a 
mapping  from  domain, 
in  the  (x,  t)  plane  to 
range,  R,  in  the  (  p  ,U) 
plane  which  belongs  to 
class  everywhere  on 
Q).  It  is  clear  that  /  is 
differentiable  at  every 
point,  such  as  c,  in 
domain,  We  note  that 
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@  and  R  are  each  in  an  space. 

Consider  the  trajectory  passing  through  point,  c,  and  another  point,  u,  in  @,  and  let  m 

be  the  distance  from  c  to  any  point,  x,  on  the  straight  line  connecting  c  and  u.  Now,  let  e  be 

the  ordered  pair  of  numbers,  (v  ,  v ),  such  that  we  =  (wv  ,  m-v )  ,  and  x  =  c+m-e . 

12  12 

Continuing,  let  D/(c)  symbolize  the  derivative  of  /  at  c,  and  note  that  D/(c)  is  a  linear 
function  at  point  c  in  Let  us  define  the  product,  D/(c)  (m  e) ,  as; 


"D 

u 

Q. 

Vj(c) 

DJicy(m-e)  =m  - 

V<.c) 

UJic) 

V,(c) 

We  note  that  existence  of  D /(c)  implies  that  the  following  inequality  must  hold: 

V  €>0,  36(e)>0,  |m-el<6(e)  =*  \fic+m-e)  -  fic)  -  iy(c)-(m-e)|  ^  e|we|  (29) 


It  follows  immediately  that: 


V  6>0,  36(6)>0,  I  we  I  <6  (e) 


— -(/(c+m-e)  -y(c)}  -  iy(c)-(e) 
m 


^  e\e 


(30) 


Inequality  (30)  proves  that  D  /(c)-(c)  is  the  directional  derivative  of  /  at  c  in  the 
direction  of  u.  It  is  not  generally  the  case  that  the  directional  derivative  of  /  in  the  direction 
of  u  is  tangent  to  the  trajectory  at  c  in  but  as  u  approaches  c  on  this  trajectory,  it  becomes 
so.  Letting  h  be  the  ordered  pair  of  numbers,  (w^,  w^,  comprising  the  components  of 

D/(c)-(c) ,  we  may  write: 


H'j(c) 

P.W 

P,(C) 

Vj(c) 

DAc)ie)  = 

Wj(c) 

— 

UJ.C) 

UJic) 

The  Jacobian  Determinant  for  the  above  transformation  should  be  written  as  seen  next. 
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(32) 


p^(c)  pfc) 
UJic)  Ufc) 


Let  us  suppose  that  J  *  0  .In  this  case,  the  following  equations  may  be  written: 


Vj(c) 

1 

t//c) 

-pf.c) 

W  (c) 

Vj(c) 

'  J 

\ 

_ 1 

P,(C) 

(33) 


It  is  clear  that  if,  and  only  if,  J  *  0  ,  Equations  (31)  and  (33)  require  that  the  linear 
function,  D /(c) ,  be  a  one-one  correspondence  mapping  the  Cartesian  space,  containing  3), 
onto  the  Cartesian  space,  containing  R. 

Let  us  introduce  that  the  following  theorem: 

Inversion  Theorem.  Suppose  function,  /  is  in  class  C^  on  a  neighborhood  of  c  e 
mapping  K'’  to  W.  Let  us  further  suppose  that  D /(c)  is  a  one-one  map  of  Rp  onto  K'’.  Then 
there  exists  V(U)  such  that  U  is  a  neighborhood  of  c,  and  V  is  a  neighborhood  of  /(c) ,  and 
/  is  a  one-one  mapping  of  U  onto  V,  and  f  has  a  continuous  inverse  function,  g,  defined  on 
set  V  to  set  U.  Moreover,  g  is  in  class  C^  on  V,  and  if  y  e  K  and  x  =  g(y)  e  U,  then  the 
linear  function,  Dg(y) ,  is  the  inverse  of  linear  function,  D /(x) . 


Note  that  the  hypothesis  that  D /(c)  is  a  one-one  map  of  onto  R^  is  satisfied  for 
p  =  2  in  ovur  system.  The  Inversion  Theorem  guarantees  the  existence  of  a  function  g 

mapping  domain  R  to  domain  3)  such  that  g  is  continuous  on  R,  and  such  that  g  =  /’^  on  R. 
We  conclude,  in  addition,  that  if  d  E  R,  and  c  =  g(d)  E  D,  then  Dg(d),  is  the  inverse  of 
D /(c) .  Consider  d  E  R,  and  let  the  trajectory  in  R,  portrayed  in  Figure  3,  pass  through 
points  d  and  z  in  R,  where  d  =  f(c)  and  z  =  Jiu)  .  By  the  Inversion  Theorem,  we  are 

assured  that  c  =  g(d)  and  u  =  g(z) .  Let  /r  be  a  particular  ordered  pair,  (Wj(c) ,  W2(c)),  and 
note  that  for  point  y  on  the  straight  line  connecting  points  d  and  z  in  R,  we  may  select  a 

positive,  real  number,  q,  such  that  q‘h  =  {q"W.^{c) ,  q’yv^ic))  ,  and  y  =  d+q'h  .  It  is  thus 
apparent  that  q  is  the  distance  from  d  to  y,  and  that  h(gd))  =  h(c)  .  Using  a  construction 
similar  to  that  applied  to  D  /(c)  in  domain  3),  the  directional  derivative  of  g  at  d  in  the 
direction  of  z  in  R  may  be  shown  to  exist,  and  be  written  as: 


Dgmh) 


Vj(c) 

1 

ufgm 

-p, ((*(<«) 

w((g((^) 

Vj(c) 

“  J 

-UJi(g(d)) 

P,(Cs(<fl) 

yv^((g(d)) 

(34) 
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As  was  the  case  for  D/(c)-(e) ,  Dg(d)ih)  becomes  tangent  to  the  trajectory  shown  in  R  as 
point  z  approaches  point  d  on  that  trajectory. 

One  may  deduce  from  Figure  3  and  Equation  (34)  that  when  z  is  located  such  that  the 
line  connecting  d  and  z  is  parallel  to  the  axis  of  abscissas,  then  =  0  ,  and  when  the  line 

connecting  d  and  z  is  parallel  to  the  axis  of  ordinates,  w  =  0  .  In  case  w  =  0  ,  Equation 
(34)  implies  that: 


Ufc) 


X  (d)  = 

(35.i) 

p 

J 

t(d)  = 

P 

U{c) 

(35.ii) 

J 

In  case  =  0  ,  Equation  (34)  implies  that: 

P/c) 

- T 

p  W 

t  (d)  = 

Equations  (35.i),  (35.ii),  (36.i),  and  (36.ii)  yield  the  relations  between  the  components  of 
D/(c) ,  and  Dg(d)  whenever  J  *  0  .It  immediately  follows  that  function  g  solves  the  pair 
of  first-order  partial  differential  equations,  seen  next,  everywhere  on  R  if,  and  only  if, 

J  ^  0  . 

X  -  U't  +  p't  =  0  (37.i) 

u  u  p 

P’x  -  p'U't  +  ,  =  0  (37.ii) 

p  p  ^ 


(36.i) 

(36.ii) 


Equations  (37.i)  and  (37.ii)  are  linear,  as  opposed  to  Equations  (27.i)  and  (27.ii),  which  are 
not.  Initially,  this  fact  would  seem  to  be  important;  but  it  turns  out  to  be  simply  a  curiosity. 

In  order  to  effectively  solve  Equations  (27.i),  (27.ii);  and  (37.i),  (37.ii),  we  must 
discuss  the  properties  of  these  equation  systems.  Consider  the  derivatives  of  /  at  point  c  in 
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and  g  at  point  d  in  R,  and  let  us  write  Equations  (27 i),  (27.ii),  and  (37.i),  (37.ii)  in 
matrix  form.  We  arrange  these  equations  in  a  style  to  accommodate  this,  as  shown  below. 
The  form  for  Equations  (27.i),  (27.ii)  is: 


+  {p^+  0}  =  0 


& 

— -p  +U‘U 

p  X  J. 


+  {0  +  U^}  =  0 


(38) 


The  corresponding  forms  for  Equations  (37 .i) 


^2  p  ^2  p 


+  {0  +  t  } 


{0.  y 


.  (37.ii)  are: 
=  0 

=  0 


(39) 


Note  that  the  partial  derivatives  of  /  at  point  c  in  ®  may  be  represented  as  the  column 
vectors  listed  following: 


— (c)  = 

9(.c) 

dx 

V^ic) 

and. 

P,(C)' 

V(f) 

(40.i) 


(40.ii) 


The  corresponding  partial  derivatives  of  g  at  point  d  in  R  are  also  represented  as  column 
vectors,  written  as: 


xM 

P 

P 


(41.i) 
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(41.ii) 


and. 


X  (d) 

t  (d) 

i/'-  '' 


Consider  point  c  in  domain  Q>,  and  point  d  in  range  R  corresponding  to  c  under  /  .  Equations 
(38)  and  (39)  motivate  the  definition  of  the  following  matrices  at  these  points: 


t/(c) 


A(c) 


Pic) 

U{c) 


(42) 


and, 

m 


0 


1 


(43) 


We  also  require: 

1  0 
I  = 

0  1 


(44) 


and, 

Nid) 


0  1 


( ll 

1  1 

- 

s 

1 

—'U 

Q. 

Ip  J 

(45) 


Finally,  we  define  matrix  ^  as: 
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(46) 


<|)  = 


0 

0 


Then  Equations  (27.i),  and  (27.ii)  may  be  written  as  the  following  matrix  equation  at 
point  c  in  domain 


=  4, 


(47) 


The  corresponding  matrix  equation  for  Equations  (37.i),  and  (37.ii)  at  point  d  in  R  are: 


+  N' — 
dU 


=  4, 


(48) 


Note  that: 

DMie)  =  |^(c)-v  (c)  +  |^(c)-v  (c)  (49) 

ox  ^  at  ^ 


and: 

Dg(d)'{h)  =  +  ^(d)’w(g(d)) 

dp  1  01/  2 


(50) 


where  v^,  are  real  numbers,  and  are  the  components  of  e  in  S).  The  real  numbers, 

are  die  corresponding  components  of  A  in  R.  Letting  M  be  a  matrix  of  elements, 

(m..)j  multiplying  the  terms  of  Equation  (47),  and  M  be  a  matrix  of  elements,  (»»»), 
multiplying  the  terms  of  Equation  (48),  we  note  that  Equation  (47)  is  equivalent  to: 
df  df 

V  ’M’A'—'v  +  V  'M'—'v  =  4>  (51) 

2  ax  1  1  0t  2  ^ 
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(52) 


and  that  Equation  (48)  is  equivalent  to: 


- dg 

w  'M’A’ — 'W 
2  dp  1 


_ dg 

+  w  ’M‘N’ — •w  =  (j) 
1  dU  2 


In  the  event  that  M  can  be  found  such  that  its  elements,  are  not  necessarily  0,  and  such 

that  M  is  not  singular,  and  v^-M‘A  =  =  B  ,  then  the  directional  derivative  of  /  in 

the  direction  of  u  from  c,  solving  Equations  (27.i),  and  (27.ii),  on  ®  satisfies: 

B'DMie)  =  4)  (53) 


If  Equation  (53)  holds,  and  M  is  not  singular,  then  the  direction  of  u  from  c  is  characteristic 
of  Equation  System  (27.i),  and  (27.ii),  and  the  differential  equation  solved  by  /  on  a 
trajectory  in  that  direction  is  ordinary.  Since  M  is  a  matrix  of  rank  2,  there  are  two  linearly 
independent  directions  in  in  which  u  can  be  located  with  respect  to  c  which  are 
characteristic. 

In  similar  fashion,  if  one  can  find  a  matrix,  M,  such  that  w^M-A  =  w^M’N  =  B  , 

and  M  is  not  singular,  then  the  directional  derivative  of  g  in  the  direction  of  z  from  d, 
solving  Equations  (37.i),  and  (37.ii),  on  R  satisfies: 

B'Dg{d)‘0i)  =  (J)  (54) 


Matrix,  M,  is  also  of  rank  2;  hence  there  are  two  linearly  independent  directions  in  which  z 
can  be  located  with  respect  to  d  in  R  which  are  characteristic,  as  well. 

The  condition,  v^M'A  =  v^M  =  B  ,  requires  that: 


/ 


m  *v  ’V  • 

^2  12  *^2 


m  *v  *U’^m  ‘V  “ 
"  21  2  22  ^2 


mii-v, 


Vx 


Vx 


Vx 


(55) 
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Corresponding  to  Equation  (54),  the  condition,  w^’M’A  =  w^M'N  =  B  ,  requires  that: 


L  V 


^  P  — ^ 

w  • — ‘tn 
2  -2  11 

V  ^ 

^  P  — ' 


1  _ 

,  w  •— -m 

I '  p 

> 

1 _ 

W  ♦ — 'TO 

1  p  22 


‘  P  _  _ '' 

-w  * — ‘U‘m  +wm 
2  ^2  11  2  12 


\ 


-w  • — ‘U‘m  +wm 
2  ^2  21  2  22 


\ 

/ 

1  — 1 

w  -m 

“  w  • 

—•U'tn 

^2 

/ 

1  11 

\ 

1 

P 

^l 

\ 

/ 

1  „ — 1 

w  -m 

-  w  - 

— 'U’m 

/ 

1  21 

1 

o. 

(56) 


Equation  (55),  in  fact,  yields  the  constraints  upon  (m^.)  if  M  is  not  to  be  singular.  First, 

Equation  (55)  requires  that  these  constraints  satisfy  Equation  pairs  (57)  and  (58),  seen  next. 
The  conditions  imposed  for  and  are: 


mii(v2-£/-Vi)+mi2-V2--  =  0 
mij-V2*p+mj2-(v2-l/-Vj)  =  0 


(57) 


The  corresponding  conditions  imposed  for  and  nijj  are: 
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(58) 


C 

m2j-V2*p+in22-(V2-t/-Vi)  =0 


Equation  (56)  yields  the  constraints  upon  {m..)  if  M  is  not  to  be  singular.  This  equation 

requires  that  the  elements  of  (ni,y)  satisfy  Equation  Pairs  (59)  and  (60)  following.  The 
conditions  on  and  are; 


—  P 

722  *W  • - 

11  2  ^2 


_  1 


m 


11 


W' — ‘U  +  w, 


\  _ / 

-  '”l2- 


1 

w,  +  W'—'U  I 
2  1  p  j 


=  0 


=  0 


(59) 


The  corresponding  conditions  on  and  Wjj  are: 

_  1 


'«2i'^2'-^ 


^21* 


"‘22'^l'- 


'  p  ^ 

1  ] 

W‘—-U  +  Wj 

-  '”22- 

w,  +  w- — 'U 

JL  ^  ^ 

^  J 

p  j 

=  0 


=  0 


(60) 


It  is  clear  that  M  is  not  constrained  to  be  singular  if  and,  only  if,  the  following  determinant 
vanishes: 
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(Vj-lZ-Vj) 


0 


(61) 


VP 


(•V2-t/-Vj  ) 


On  the  other  hand,  matrix  M  is  not  constrained  to  be  singular  if  and,  only  if: 


^  p  ' 

w  •— -17  +  w 

1  1 

C 

\  / 

P  J 

0 


(62) 


Note  that  Equation  (61)  requires  that: 


v^'(U±C)  -  Vj  =  0  (63) 

Suppose  V2*(  U  +  C)  -Vj  =  0.1n  this  case,  (v^'U-v^)  =  -v^'C  ,  and  Equations  (57) 
require: 


m  =  m  • — 

12  11  c 


(64) 


Matrix  M  will  be  singular  unless  we  incorporate  the  other  form  of  Equation  (63)  in  Equations 
(58).  The  other  form  is:  V2'(  U-C)  -Vj  =  0,  which  means  that  (y^U-v^  =  v^C  ; 
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whence: 


m  = 


21 


(65) 


Let  =  Wjj-Vj  and 


B  may  now  be  expressed  as: 


(66) 


It  is  apparent  from  the  form  of  matrix  B  that  its  elements  are  associated  with  the 
characteristic  directions  in  S*  as  illustrated  Table  1. 

Table  1 


fc:  V2-(l/  +  C)  -  Vj  =  0 

b^:  v^iU-C)-v^=0 

In  each  of  the  characteristic  directions  associated  with  b^  and  b^  respectively,  Matrix 

Equation  (53)  yields  a  linearly  independent  form  depending  upon  D/(c)*(c) ,  written  as 
follows: 


/  \ 

/ 

•V 

+  b- 

p 

+  — -t/ 

•V 

C  ^ 

1 

I 

C  f 

2 

\ 

\ 

J 

(  C 

N 

c 

\ 

.  -  — *0  + 

u 

•V  + 

b  • 

— *p. 

+  u 

i  p  * 

X 

1 

II 

p  ' 

t 

/ 

=  0 


*v  =  0 
2 


(67) 


On  the  other  hand,  Equation  (62)  requires  that: 
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(C'Wj  -  p'Wj)  •  (C‘Wj  +  P'Wj)  =  0 


(68) 


Equation  (68)  implies  either  (Ow^-p-Wj)  =  0  or  (C'Wj  +  p'Wj)  =  0.  In  the  case  of 
matrix  M,  if  the  first  condition  holds,  then: 


m 


12 


(70) 


As  in  the  case  of  matrix  M,  M  will  be  singular  unless  the  second  condition  is  applied  to 
define  the  elements  of  its  second  row.  Consequently, 


(71) 


Using  Equations  (70)  and  (71),  we  may  write  Equation  (56)  as: 


w  -M'A  =  w  'M'N 
2  1 


u 


+  m 


11 


JP 

C 


*w 


f  —  0 

_  u  _ 

-m  •  — 

•w 

m  ' —  +  m 

22  C 

2 

22  C  22 

‘W 


(72) 


Defining 


and 


matrix  B  may  be  written  as: 
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B  = 


F.i 

n  C 


-b 


-b 


U 


-1 


bj'iU  +  C) 


(73) 


Hence,  the  characteristic  directions  for  g  =  /  \  solving  Equations  (37.i)  and  (37.ii),  may  be 

associated  with  b  and  b  as  follows  in  Table  2: 

/  n 


Table  2 


b  : 

(C'Wj  +  p'Wj)  =0 

/ 

(C-Wj  -  p'Wj)  =0 

In  this  case,  Matrix  Equation  (54)  yields  a  linearly  independent  form  in  each  of  the 
characteristic  directions,  depending  upon  Dg(d)-(h),  which  may  be  written  as: 


bj[-Xp+(,U+C)'tp]-w^  +  b^’[-Xfj  +  (U+C)-tu]-w^  =  0  (74) 

'^'[-Xp+{U-C)‘tplw^  +  ^‘[-Xu  +  (U-C)'tuyw^  =  0 


Now  consider  a  particular  point,  p  ,m  S)  and  let  p  be  its  image  under  /  in  R. 
Without  loss  of  generality,  we  may  select  the  direction  in  S  associated  with  b^.  Since  f 

belongs  to  class  C\  and  p  is  not  special,  one  may  apply  Inequality  (30)  to  D /(p)-(e)  solving 
Equation  (53)  to  construct  a  unique  curve,  ,  passing  through  point  in  the  selected 
direction.  The  condition  that  /  belongs  to  class  C‘  implies  that  Inequality  (30)  may  be  used 
to  construct  F^  as  the  limit  of  a  Cauchy  sequence  of  continuous  curves  passing  through 
point  p  .  The  limit  is  the  unique  curve,  F^ ,  passing  through  point  p  ,  which  has  a 
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continuously  turning  tangent  in  S.  Curve  is  a  characteristic  curve  of  Equation  System 
(27.i)  and  (27.ii).  It  is  apparent  that  Equation  (53)  admits  the  construction  of  precisely  one 
other  curve,  Tjj,  passing  through  p  ,  which  is  also  a  characteristic  curve  of  Equation  System 

(27.i)  and  (27.ii).  Let  us  associate  F^  with  b^,  and  F^  with  b^,  in  matrix  B  in  the 
subsequent  discussion. 

Characteristic  curves,  F^  and  Fjy,  have  images,  T/  in  R  under  /  ,  which 

intersect  at  p  .  Both  of  these  pairs  of  curves  are  illustrated  in  Figure  4.  Here,  and  ajj 
are  linearly  independent  tangent  vectors  associated  respectively  with  F^  and  F^.  Clearly, 
curves  and  are  characteristic  curves  of  Equation  System  (37.i)  and  (37 .ii)  provided 

that  J  0  .  In  this 
case,  a  I  and  ttu 
are  linearly  in¬ 
dependent  tangent  ^  U 

vectors  associated 

respectively  with 

and  Yjj-  Obviously,  F^ 

and  Fj^  are  the  images 

of  T/  F//  ® 

under  g  ,  when  g 
exists. 

Note  that  the 

existence  of  g  =  /  \ 
te.  J  *  0  ,  is  not 
essential  to  the 
construction  of  F^  and 

F^  in  This 

development  leads  to  a 
class  of  special  wave 

functions  which  solve  Equation  (53),  but  do  not  require  the  condition,  J  *  0  .  The 
expansion  of  a  perfect  gas  from  an  infinite  reservoir  into  a  void  main  is  described  by  one  of 
this  class  of  special  wave  functions,  which  is  described  in  Section  5. 

Continuing  with  the  condition  that  J  *  d  ,  let  and  O ^  be  distance  parameters 

associated  with  F^  and  F^,  respectively;  and  let  Oj  and  Ojj  be  distance  parameters 

associated  with  T/  F//-  In  domain  3>,  we  may  make  the  following  table  of  associations 

between  the  characteristic  directions  and  derivatives  with  respect  to  the  distance  parameters: 


25 


Table  of  Associated  Characteristic  Directions  and  Total  Derivatives 


r  : 

I 


and 


dx 

la 

I 

dt 

^2  ~  H 

I 

dx 

V  =  - 

1  da 

n 

dt 

V  =  — — 

2 


Corresponding  Associations  in  Region  R 


r : 

I 

r/ 


and 


w 

1 


dp 

do 

n 


w 


dU 

d^ 

u 
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Now  at  point  p  in  domain  S>,  Table  1,  found  on  page  22,  demonstrates  that  the  domain 
points  of  function  /  solving  Equations  (27 .i)  and  (27.ii)  satisfy  the  ordinary  first-order 
differential  equations  shown  in  Table  3,  in  the  indicated  characteristic  directions. 

Table  3 


(t/  +  C)-—  =  0 
da 

/ 


Equations  (67)  yield  first-order,  ordinary  differential  equations  for  the  corresponding  range 
points  of  /  in  the  directions  of  and  Fj^.  These  Equations  are  illustrated  in  Table  4. 

Table  4 

F  :  iL  +  =  0 

^  da  C  da 

I  I 

F  :  _  P .  ^  Q 

"  do„  c’do„ 

n  n 

In  case  J  0  ,  function  g  exists  such  that  g  =  f  ^  •  Set  R,  in  the  (  p  ,U)  plane,  is  then  in 
the  domain  of  g  ,  and  g  solves  Equations  (37.i)  and  (37.ii)  on  R.  At  point  p  in  set  R, 
Table  2,  found  on  page  24,  demonstrates  that  the  domain  points  of  function  g  solving 
Equations  (37. i)  and  (37.ii)  satisfy  the  ordinary,  first-order  differential  equations  shown  in 

Table  5,  below,  in  the  characteristic  directions,  and  R- 


Table  5 


p  dU 
C  da 


p  .  “P  P  « 

da  C  da 
n  n 

Equations  (74)  yield  first-order,  ordinary  differential  equations  for  the  corresponding  range 
points  of  g  in  the  directions  of  and  p^.  These  Equations  appear  in  Table  6  following. 
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Table  6 


dx 

do 

I 


dt 

(t/  +  C)-—  =  0 
do 

I 


dx  dt 

-  {u-cy-^  =  0 

n  n 


The  Second  Law  of  Thermodynamics  and  the  Perfect  Equation  of  State  for  a 
calorically  perfect  gas  may  be  embodied  in  Tables  4  and  5  to  produce  particularly  simple 
first-order  differential  equations.  In  each  characteristic  direction  in  at  a  point,  we  write  the 
following  associated  forms,  noting  that^  R  is  the  gas  constant  for  the  Perfect  Equation  of 
State  in  Tables  7  and  8: 

Table  7 


piOj)  =  piOjyR-TiOj) 


ds  _  y  1  dp 

dOj  P  ^^1 


y  ^  dT  1  dp 

Y-1  dOjj  p  dOjj 


In  each  characteristic  direction  in  set*  R,  the  domain  of  function 
following  corresponding  forms: 


Table  8 


g 


,  at  a  point,  we  write  the 


r^:  =  piOjyR’TiOj) 


r  •  p(°ii>  =  piOjjyR  TiOjj) 

//■ 


ds  y  dT  1  dp 

T  —  =  —'R  —  - 

dOj  Y“1  dOj  P  dOj 


ds  y  dT  1  dp 

r =  —  i?*—  -  — •-= 

do^  Y-1  do„  P  do„ 


^  R  refers  to  the  range  of  /  except  when  used  in  the  context  28  of  the  Perfect  Equation  of  State. 


From  om  hypothesis  that  the  Euler  Equations  of  Motion  hold,  it  has  been  shown  that  the 
wave  is  isentropic.  Hence,  the  thermodynamic  relations  of  Table  7  result  in  the  following 
ordinary,  first-order  differential  equations  for  the  range  points  of  /  ,  holding  in  each  of  the 
characteristic  directions  in  &  : 


Table  9 


r  : 

dp 

2 

p  dC 

dp  p  dU 

/ 

_  _ 

Y-1 

C  dOj 

da^  C  da^ 

r„: 

dp 

2 

p  dC 

dp  p  dU 

// 

~^n  ' 

Y-1 

C 

da  C  da 

n  n 

For  completeness.  Table  10  illustrates  the  ordinary,  first-order  differential  equations  in  each 
characteristic  direction  in  set  R,  the  domain  of  inverse  function,  g  . 


Table  10 


dp 

2 

p  dC 

do 

p  dU 

r : 

/ 

Y-1 

— =  + 

da 

I 

C  da^ 

F : 

2 

p  dC 

dp 

p  dU 

II 

Y-1 

C  dTj^ 

da~ 

n 

C  d^ 

tt 

Through  Tables  1  -  10,  one  can  see  that  function  /  ,  solving  Equations  (27  i)  and 
(27.ii),  is  transformed  through  a  one-to-one  correspondence  to  characterize  the  solutions  to  the 
following  set  of  ordinary,  first-order  differential  equations  which  appear  in  Table  11: 
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Characteristic 


Table  11 
Domain  @ 


Range 


dx 

1 

dt 

-  (1/  +  C)*  -  0 

do 

/ 

2 

Y-1 

dC 

.a/ 

dU 

daj 

dx 

dt 

2 

dC 

dU 

-  {u-cy—  =  0 

do 

n 

Y-1 

da^j 

The  corresponding  ordinary,  first-order  differential  equations  applying  to  inverse  function  g 
are  illustrated  in  Table  12;  again,  for  the  sake  of  completeness. 


Characteristic 

f : 

/ 


Table  12 
Range  9 


dx 

da 

I 


dt 

-  (1/  +  C)*—  =  0 
da 

I 


dx 

Jo” 

n 


dt 

-  =  0 

“’’n 


Domain 

2  dC  dU  ^ 

- +  __  =  0 

Y“1  daj  daj 

2  dC  dU 

- -  __  =  0 

Y“1  da„  da„ 

u  u 


Letting^  a(a^)  ,  P(op  ;  and  a(o^)  ,  p(o^^)  ,  be  strict  functions  of  the 

distance,  or  curve  parameters  of  characteristic  curves,  ,  and  F^  ,  and  requiring  that 
each  of  these  distinct  function  pairs  be  defined  at  each  point  on  each  of  the  characteristic 


^  R  is  range  of  f  transformed  according  to  the  relations  of  Table  7. 
^  a(o  )  *  a(o  )  »  P(<y )  ’ 


curves,  ,  and  F^  ,  indefinite  integrals  may  be  written  on  these  curves  in  the  form 


seen  in  Table  13: 

Tabie  13 

Characteristic 

Domain  @ 

Range  R* 

F  : 

7 

X  -  J(u  +  cydt  = 

“(V 

J---C  *  V  =  P(o„) 
y-l 

F  : 

77 

X  -  j(U-C)'dt  = 

«(o,) 

-i--C  -  V  =  P(o,) 

Y-l 

Table  13  may  be  thought  of  as  the  restriction  of  /  to  F^,  written  as  /  I"  F^  ,  in  the  first 

row,  and  as  the  restriction  of  /  to  F^,  written  as  /  f  Tjj  ,  in  the  second  row. 

The  corresponding  result  for  g  =  f  ^  is  shown  in  Table  14,  where^  Ct(o])  » 

;  and  a(^)  »  P(^)  >  corresponding  strict  functions  of  the  distance,  or 

curve  parameters  of  characteristic  curves,  ,  and  ^ 


Table  14 


Characteristic 

Range 

Domain 

X  -  J(u  +  c)‘dt  =  a(^) 

2 

^_J-C  .  U  =  £(a„) 

^77- 

X  -  j(u-c)’dt  =  a(^) 

^•C  -  1/  =  P(o,) 

*  R  is  the  range  of  /  transformed  according  to  the  relations  of  Table  7. 


a(op  >  a(op 


P(o")  »  P(5~)  respectively,  distinct  pairs  of  functions 
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Table  14  may  be  thought  of  as  the  restriction  of  g  to  P/  >  written  as  g  I  P/  >  ^^st 

row,  and  as  the  restriction  of  g  to  p^,  written  as  g  1  p^^  ,  in  the  second  row. 

Noting  that  for 
each  point  p  in  &,  the 

domain  of  /  solving  ^  u 

Equations  (27. i)  and 
(27.ii),  has  two 
characteristic  curves, 

P^  and  passing 

through  it  in  linearly 
independent  directions, 
and  a  single  pair  of 
distance  or  curve 
parameters,  Oj  and  Ojj, 

associated  with  it,  one 
may  create  a  mapping 
from  the  (x,  t)  pair  to 
the  ( Oj,  Ojj)  pair 
associated  with  each 
point  in  @.  Let  9^  be 

thesetof(aj,  Figures 

pairs  under  this 

mapping.  This  is  an  injective  mapping,  written  symbolically  as  : 

S!  ’4  ^  (75) 

Referring  to  Figure  5,  suppose  we  choose  to  move  along  the  characteristic  curve  passing 
through  p  .  In  this  case,  cmve  parameter  Gj  varies  with  position  on  F^ ,  but  G remains 
unaffected,  since  the  other  characteristic  direction  is  linearly  independent. 

Let  =  jUW^,  and  let  be  the  range  of  J\f .  Then  Dmn  D =  Fj  and 

Bng  t/f^D  ^  .  Using  set  builder  notation,  t/f'  may  be  described  as: 

^  =  {(n,  m)  \  [n  =  (x,  t)  A  n  e  Tj  A  m  =  (Oj,  o^)A  m  e  ]  A 
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Clearly,  must  have  the  property  that  V  n^,  ^  ^ 

not  to  be  simply  a  fixed  point.  Evidently,  cc(a^)  P(a^j)  have  fixed  values  on 

characteristic  Fj .  We  could,  as  well,  have  chosen  characteristic  curve  to  construct  ; 
hence,  we  conclude  that  a(Oy)  and  P(ay)  have  fixed  values  on  characteristic  Fyy,  also. 

Turning  again  to  Figure  5,  let  us  consider  point  p  in  R,  the  image  of  point  p  in  @ 
under  /  .  Now  R  is  the  domain  of  g  ,  and  so  p  has  two  characteristic  curves,  and  > 

passing  through  it  in  linearly  independent  directions,  and  a  distinct  pair  of  distance  or  curve 
parameters,  ^  and  associated  with  it,  just  as  was  the  case  with  p  in  0.  Choosing  Ty  > 

and  letting  be  the  set  of  pairs,  (Oy,  Oyy),  associated  with  the  pairs,  ( p ,  17 ),  along  Py^ 
we  may  write  the  following  injective  mapping  on  R: 

^ 

Proceeding  in  the  same  manner  with  function  g  as  we  did  with  function  /  ,  let 
=  rA(\F^,  and  let  ^  range  of  J\f .  Then  Dmn  R  =  Fj  and 

RngJ^  n  '^  =  .  Using  set  builder  notation,  may  be  described  as: 

(n,  w)  1  [n  =  (p,  t7)  A  n  e  r|  A  m  =  (^,  A  m  e  ^  ]A[^=7^  11 

We  find  that  has  the  property  that  V  «  ,  n  ,  it ,  n  e  p  =*  (O/?),  =  if  F„  is 

not  to  be  a  fixed  point.  This  shows  that  a(a^)  and  P(^)  have  fixed  values  on 
characteristic  Obviously,  we  have,  without  loss  of  generality,  shown  that  a(o^)  and 
P(^)  have  fixed  values  on  characteristic  p^,  in  addition.  We  note  here  that  P(ay)  and 
sometimes  called  the  Riemann  Invariants. 

Continuing,  we  shall  refer  to  Fy  and  Fy  as  characteristics  of  Type  I  in  and  we 

shall  refer  to  Fyy  and  Fyy  as  characteristics  of  Type  II  in  R.  The  pairs  of  Type  I  and  II 

characteristic  curves  form  an  untangled  net  covering  domain  ®  of  /  ,  as  shown  m  Figure  5. 
The  pairs  of  Type  I  and  II  characteristic  curves  form  an  untangled  net  covering  domain  R  of 

g  =  f  ^ ,  when  it  exists.  This  corresponding  net  is  also  shown  La  Figure  5. 
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5. 


SIMPLE,  ONE-DIMENSIONAL  WAVES 


Our  characterization  of  simple  waves  begins  with  the  following  definition. 


Definition.  Suppose  the  there  exists  a  domain  S)  in  the  (x,  t)  plane  for  which  /  solving 
Equations  (27.i)  and  (27 .ii)  belongs  to  class  C*.  *  Let  us  suppose  that  all  of  the  characteristics 
of  a  single  type  covering  e.g.  Type  I,  are  straight  lines.  Then  if  R  is  the  plane  of  ( p ,  U) 
pairs,  we  say  that  this  function,  symbolized  by: 


is  a  simple  wave. 


Theorem.  A  function  /  which  belongs  to  class  and  which  solves  Equations  (27 .i)  and 
(27.ii)  on  domain  in  the  (x,  t)  plane,  is  a  simple  wave  on  &  if,  and  only  if.  Determinant  J, 
defined  by  Equation  (32),  vanishes  at  every  point  in  3). 

Proof.  Consider  the  directional  derivative  of  /  at  point  p  in  3,  and  recall  its  form  from 
Equation  (31)  in  Section  4.  We  write: 


(77) 


Also  recall  that  the  determinant  of  D  /(c)*(e)  is: 


P  P 
U  U 


(78) 


^  A  function  on  a  domain  Q)  in  is  of  class  C*  if  all  of  its  partial  derivatives  exist  and  are  continuous  on  0). 
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Now,  by  Cramer’s  Rules,  the  following  Equations  must  hold: 


w  p 
w  U 

2  t 


P 

X  1 

U 

X  2 


(79.i) 


(79.ii) 


Suppose  J  =  0  ;  but  let  us  require  that  an  ordered  pair,  (Vj,  Vj)  ,  exists  at  every 
point  in  In  this  case,  the  following  identity  holds,  provided  that  (o  exists  : 


Pr 


(d 


(80) 


Now  G)  exists  everywhere  in  2)  if,  and  only  if,  0  at  every  point  in  2);  which  will  be 

the  case  if  U  depends  upon  p  everywhere  in  2.  But  we  are  assured  by  Equation  (24.iii) 
that  this  condition  holds  .  Also,  since  Equations  (27 .i)  and  (27.ii)  hold  everywhere  in  2,  the 
condition  that  J  =  0  ,  when  applied  to  them,  guarantees  that: 

CD  =  ±  -  (81) 

P 


Without  loss  of  generality,  we  may  set  CD 


+  —  at  point  p  in  2.  Continuing,  the 

P 

condition  that  J  =  0  ,  when  applied  to  Equations  (79.i)  and  (79.ii),  requires  that: 


w  'U  -  w  -p  =0  (82.i) 

w  'U  -  w  'p  =0  (82.U) 

\  t  2*^t 
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We  reiterate  that  the  isentropic  condition,  guaranteed  by  Equation  (24.iii),  implies  that  (O 
exists  everywhere  in  Hence,  Equations  (82.i)  and  (82.ii)  reduce  to  the  single  condition: 


w 

1 


—  -w  =  0 
C  2 


(83) 


Let  a  be  a  general  distance  parameter  in  the  direction  specified  by  c  in  D/(c)-(e) , 
and  let  us  associate  the  following  derivatives  with  (Vj,  Vj)  and  (Wj,  w^)  at  a  point  in 
domain 


dx 

_  dp 
da 

do 

dt 

dU 

w,  =  — 

do 

^  da 

Then  J  =  0  implies  that  Equation  (83)  has  the  following  invariant  form,  for  any  direction 
specified  by  e  in  D /(c)-(c) ,  at  point  />  in 


dp  p  dU 
da  C  da 


(84) 


We  could,  as  well,  have  made  the  choice^  (o  =  -  —  . 

P 

Recall  from  Section  4  that  the  existence  of  g  =  f  ^ ,  i.e.  J  *  0  ,  is  not  essential  to 
the  construction  of  and  F^  in  Hence,  we  may  align  Df(c)ie)  at  point  in  with 

the  characteristic  solutions  of  Equation  (53),  which  can  be  constructed  even  when  J  =  0  . 
When  we  do  this,  the  following  associations  appear  in  the  characteristic  directions,  as 
illustrated  in  Table  15  : 


•  The  form  of  Equation  (84)  would  be 


da 


9  dU  . 

—  • —  =  0  in  this  case. 
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Table  15^ 

Characteristic  Associations  When  J  =  0 


dx 

da 

I 


(u  +  cy 


=  0 


dx 


dt 

-  (u-c)‘—  =  0 

da 

n 


dp 

p  dU 

=  0  * 

dp  _ 

p. 

dU 

=  0  * 

daj 

~claj 

dajj 

c 

dajj 

dp 

p  dU 

=  0 

dp 

p. 

dU 

=  0 

la 

/ 

~C  la 

I 

do 

u 

c 

da 

n 

Both  ordinary,  first-order  differential  equations  for  the  range  points  of  f  on  apply,  so 

we  conclude  that  on  this  Type  I  characteristic,  the  range  points  of  /  satisfy  the  conditions 
shown  below: 

iP  =  0  P 

U  =  Hiaj^) 


where  GiPjj)  and  HiOjj)  are  strict  functions  of  the  curve  parameter  for  a  characteristic 
of  Type  II  intersecting  F^  . 

'  The  *  appearing  in  Table  15  indicates  the  form  of  Equation  (84)  when  D  /(c)‘(c)  is  aligned  with  the  indicated 
characteristic,  and  J  =  0  ,  37 


Since  the  characteristic  directions  at  point  p  are  linearly  independent,  curve 
parameter  is  fixed  at  every  point  on  .  Since  point  p  was  chosen  arbitrarily,  it 

follows  that  when  J  =  0  ,  G(o^)  and  H(o^)  have  fixed  values  all  along  passing 

through  point  p  .  In  view  of  the  condition  that  the  wave  is  isentropic,  we  refer  to  Table  7 
and  find,  consequently,  that  —  =  0  .  It  follows  that  C  =  JvWT  is  necessarily  a  strict 

dOj  ^ 

function  of  on  F^  .  Hence,  Table  15  shows  that  characteristic  F^  is  a  straight  line. 

Letting  aiid  B(o^)  be  two  more  strict  functions  of  a  curve  parameter  for 

a  characteristic  of  Type  II  intersecting  F^  ,  the  equation  for  F^  ,  and  the  equation  for  the 
Riemann  Invariant  on  F^  appear  on  the  first  line  of  Table  16.  The  corresponding  forms  for 
the  Type  I  characteristic  appear  on  the  second  line: 

Table  16 

Characteristic  Domain  3)  Range  R 

-L-c  +  V  =  B(V 

Y_1  u 

^•C  -  U  =  BiOj) 
Y-1 


X  -  (l7  +  C)-r  =A{ajj) 


n 


X  -  j‘(U-C)'dt  =  A(aj) 


Clearly,  A(aj)  and  ;  and  AiOjj)  and  B(o^)  have  fixed  values,  as  well,  on  the 

indicated  characteristics',  in  their  associated  rows,  in  Table  16. 

On  the  other  hand,  suppose  F^  is  a  straight  line  in  S*.  It  can  be  shown,  as  a  result, 

that  =  0  and  =  0  everywhere  on  F  .  Consider  a  particular  point  p  on  F  , 

dOj  dOj  '  ‘ 

and  note  that  the  following  matrix  equation  must  hold: 

Aia^)  ,  A(o^)  ;and  B(op  ,  B(o^  )  are,  respectively,  distinct  pairs  of  functions. 


1 


(85) 


Ufp)  -plip) 

w  (p) 

Vj(p) 

. 

=  J- 

-U(p)  pip) 

V,0>) 

where  the  pairs  ,  (v^,  v^)  and  (w^,  w^)  ,  are  defined  on  page  12.  By  hypothesis, 


passes  through  p  in  S);  thus  one  of  the  pair,  (Vj,  v^)  ,  must  be  different  from  zero.  If  we 

dp  dU 

require  J  *  0  ,  then  Equation  (85)  implies  that  w  =  — —  and  w  =  -y-  ;  but  either 

^  dOj  dOj 

|wjl  >  0  or  |Wjl  >  0  ,  in  this  case,  contradicting  our  supposition.  Hence  J  =  0  ,  and 


our  theorem  is  proved.  ■ 
Each  of  the 
characteristics  of  Type 
I  form  a  straight  line 

for  our  choice  of  (0  in 
Equation  (81).  A  sketch 
of  the  characteristic  net 
for  such  a  wave  is 
shown  in  Figure  6.  The 
range  of  /  on  Q)  is 
the  single  curve  in  the 
( p ,  U)  plane,  depict¬ 
ed  as  R  in  Figure  6. 
Every  point  on  the 

illustrated,  such 

as  points  p^,  p^,  p^, 
and  p^ ,  maps  to  the 


U 


single  point,  p  ,  on 
curve  R,  for  the 

example  shown  here.  Figure  6 

A  simple  wave 

whose  straight-line  characteristics  in  domain  3!  are  of  Type  I  is  a  forward-facing  wave.  If  we 


had  chosen  (O  =  -  —  at  point  p  in  the  straight-line  characteristics  would  have  been 

P 

of  Type  II,  and  the  simple  wave  would  have  been  backward-facing.  If  all  of  the  straight-line 
characteristics  intersect  at  a  point  on  the  bound  of  the  wave  is  centered. 

The  head  of  the  wave  of  interest  in  this  analysis  is  located  at  plane  0-0  in  Figure  1.  It 
is  of  characteristic  Type  II.  The  flow  enters  the  wave  here  at  critical  velocity  and  proceeds 
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into  the  rarefied  region  at  supersonic  speed.  It  is  clear  that  the  thermodynamic  state  on  the 
head  characteristic  is  constant,  and  that  characteristic  is,  therefore  straight.  Obviously,  every 
characteristic  of  Type  I  in  the  domain  S)  of  this  wave  has  an  intersection  with  the  head 
characteristic  at  plane  0-0,  and  each  of  these  intersections  is  at  the  same  thermodynamic  state. 
It  is  easy  to  show  that  J  =  0  at  every  point  p  in  2),  as  a  consequence.  Hence,  the  subject 
wave  is  a  simple,  backward-facing,  centered  expansion  wave  -  sometimes  known  as  a 
centered  rarefaction  wave. 

There  is  a  shock  at  the  tail  of 
the  backward-facing  expansion  wave 
which  is  not  in  the  domain  of  the 
wave,  and  to  which  characteristic 
theory  does  not  apply.  We  need  to 
establish  that  this  shock  exists.  Since 
every  characteristic  of  Type  n  is 
straight  in  domain  one  can  show 
that  the  infimum  or  greatest  lower 
bound  of  characteristics  ,  is  also  a 
straight  line,  i.e.,  it  has  a  constant 
velocity  of  propagation.  Consider 
Figure  7  and  let  a  fibced  rectangular 

control  volume  4-  be  such  that  face  L 

of  is  to  left  of  and  such  that  face 

R  of  4  is  to  its  right.  Then  control 

volume  4  must  proceed  at  constant  velocity,  ,  which  is  the  velocity  of  propagation  of  A 

The  right  face  of  control  volume  4  is  in  vacuo.  This  being  the  case,  it  must  be  that  p  =  0 
at  R  in  Figure  7,  where  p  s5mibolizes  static  pressure  in  this  portion  of  the  discussion. 

Clearly,  L  is  interior  to  domain  hence  p  >  Q  at  face  L. 

Recall  that  the  flow  in  domain  @  of  the  expansion  wave  is  isentropic.  The  fact  that  the 
static  pressure,  p  ,  at  a  point  c  in  the  interior  of  @  is  such  that  p{c)  >  0  implies  that 

both  p(c)  >  0  and  K(c)  >0  .  If  the  gas  is  calorically  perfect,  the  thermodynamic 
properties  at  point  c  are  related  to  those  at  face  L  by  the  following  two  equations  : 


4 - S*'  Kv 


Figure  7 


(86) 


and 
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/ 


1 


(87) 


Continuing,  let  us  establish  the  following  definitions  with  reference  to  where  subscript  L 
refers  to  the  variable  value  at  face  L  in  Figure  7: 

Definition  Face  L  Face  R 

2 

Total  Energy  ^  ^  ^  ^ 

~  ^  2 

Total  Enthalpy  h  =  e  +  —  ^  ^ 

L  L 


Impulse  h^-Pi-  Pr  ® 

Mass  Flux  =  -  Pi  (^4-^cv) 

Clearly,  every  nested  sequence  of  control  volumes  {<i}i,i,oo»  by  definition  of  4^,  contains 
and  has  as  its  limit,  ^  with  faces  L  and  R  on  its  left  and  right,  respectively.  Let  us  refer  to 
this  limit  as  4^  •  We  now  observe  that  if  the  gas  is  calorically  perfect,  the  specific  internal 
energy,  and  the  sound  speed  C,  are  related  according  to  Equations  (88)  and  (89),  shown  next: 

Pl  (88) 

—  =  (Y-l)-«i 
Pl 


C  = 


ypi 


(89) 
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Now  in  the  case  of  our  backward-facing  expansion  wave,  one  can  show  that  Riemann 
Invariant  ,  from  Table  16,  is  fixed  everywhere  in  S>.  Furthermore,  for  the  limiting 

containing  i.e.,  4^^ ,  Equation  (86)  implies  that  =  0  .It  follows  from  Equations  (88) 
and  (89)  that  t/^  =  for  the  limiting  case.  This  velocity  is  the  upper  limit,  or  escape 

velocity  for  fluid  elements  in  It  is  clear  enough,  then,  that  as  well.  It  now 

follows  that  with  respect  to  4^  : 


o 

11 

( 

II 

(90.i) 

L  R 

II 

II 

O 

(90.ii) 

II 

II 

o 

(90.iii) 

(90.iiii) 

=  0 

Equations  (90.ii),  (90.iii),  and  (90.iiii)  are,  respectively,  conservation  of  energy, 
momentum,  and  mass  flow  in  the  reference  frame  of  </.  They  constitute  the  Rankine-Hugoniot 

conditions,  and  are  satisfied  with  respect  to  the  limit  control  volume,  4  m  •  Note  that  they  are 
trivially  satisfied  under  the  condition  of  no  mass  flux  at  either  face  of  the  wave,  as  would  be 
the  case  if  ^were  a  contact  surface.  This  condition  would  normally  require  that  -  Uj^  , 
making  the  convective  gas  velocity  continuous  at  the  location  of  However, 

=  0  for  4  m  containing  so  =  ntj^  =  0  even  when  ^i*  •  Since 

=  B{Ojj)  at  face  L,  and  £7^  =  0  at  face  R,  it  is  clear  that  c  4m  is  the  site  of  a 
velocity  discontinuity  -  an  abrupt  velocity  jump  from  t/  =  0  to  U  =  BiOjj)  propagating 

at  escape  velocity  for  the  expansion  wave  on  domain  Hence  is  a  shock  with  its  front 
facing  the  vacuum,  and  its  back  facing  the  expansion  wave. 

6.  CALCULATION  OF  THE  ONE-DIMENSIONAL,  BACKWARD-FACING 

EXPANSION  WAVE 


Consider  Figure  1  from  Section  3.  We  are  assuming  that  the  flow  is  one-dimensional 
everywhere,  and  that  the  reservoir  is  infinite  in  extent.  Hence,  when  the  rarefaction  is  initiated 
at  time  zero,  there  can  be  no  influence  on  the  thermodynamic  state  of  the  reservoir.  There  can 
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be  no  local  flow  effects  in  the  gas  to  the  left  of  plane  0-0  without  obviating  the  one- 
dimensional  flow  field  hypothesis.  This  is  to  say  that  the  reservoir  has  sufficient  depth  to 
maintain  a  fixed  stagnation  pressure  and  temperature  everywhere  to  the  left  of  plane  0-0. 
Consequently,  the  reservoir  is  able  to  maintain  a  constant,  perpetually  choked  condition  at 
plane  0-0  after  the  initiation  of  the  rarefaction.  Letting  be  the  stagnation  temperature  at 

plane  0-0,  and  be  the  stagnation  pressure  there,  these  thermodynamic  quantities  are 

fixed  at  plane  0-0  for  all  time  following  the  initiation  of  the  rarefaction.  The  gas  velocity  at 
plane  0-0  is  fixed  after  time  zero,  as  a  result,  maintaining  a  value  expressed  as: 


U  = 

cr 


2yRT^ 

^J  (Y+1) 


(91) 


The  head  of  the  wave  is  adjoined  to  plane  0-0,  which  is  maintained  at  a  constant 
thermodynamic  state.  Hence,  J  =  0  everywhere  in  domain  S*  of  the  wave,  with  the 
consequence  that  it  is  a  simple  wave.  Furthermore,  we  deduce  that  the  lead  characteristic, 
which  is  at  the  head  of  the  wave,  must  permanently  coincide  with  plane  0-0  after  time  zero, 
since  the  wave  cannot  penetrate  the  reservoir  and  disturb  its  thermodynamic  state.  This  means 
that  Type  II  characteristics  must  be  straight,  making  the  wave  backward-facing.  This  means 

also  that  there  exist  curve  parameters  and  for  characteristic  families  Tj  and 
;  and  functions  and  AjiOj)  ,  and  fimctions  BfOjj)  and  BJjOj)  ,  such 

that  AljOj^  maintains  a  fixed,  but  distinct  value  for  each  characteristic  Tj  ,  and 

maintains  a  fixed  value  for  all  characteristics  ;  and  AJ^Oj)  and  BjfPj)  each 

maintain  fixed,  but  distinct  values  for  each  characteristic  F^  .  In  this  case,  we  have  the 
following  correspondences  on  the  characteristics  of  the  expansion  wave  in  domain 

Table  17 

Riemann  Invariant 
^•C  .  V  = 


Characteristic  Equation  in  Domain  2) 

r:  X  -  jiu^cydt  =Afa„) 

-  (.u-cyt  =  A 
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It  is  obvious  that  the  Type  I  characteristics  span  the  wave;  hence  is 

determinate  at  plane  0-0,  the  plane  of  choked  flow,  as  well  as  having  a  fixed  value 
everywhere  in  Note  that  at  plane  0-0,  we  may  write: 

c  (92) 


In  view  of  Equation  (92)  and  the  definition  of  ,  we  find  that: 


2yiy+iyRT^ 

'i  (Y-i)" 


(93) 


Let  e  be  the  total  energy  of  the  flow  at  a  point  m  in  domain  @  and  let  p-v  represent 
the  flow  work  required  to  compress  the  flow  adiabatically  from  its  thermodynamic  state  at 
limit  velocity  to  its  actual  thermodynamic  state  at  point  m  .  Also,  let  represent  the 

flow  work  to  adiabatically  stagnate  the  flow  from  that  same  thermodynamic  state  at  limit 
velocity.  Continuing,  let  represent  the  local  stagnation  internal  energy,  and  note  that 

.  Now  the  First  Law  of  Thermodynamics  requires  the  change  in  total  energy 

between  these  two  states,  at  point  m  ,  in  our  adiabatic  system,  to  be  a  consequence  of  the 

difference  in  flow  work  required  to  reach  each  of  them  from  their  common  thermodynamic 
state  at  limit  velocity.  This  condition  is  expressed  by: 

(e^  -  g)  +  (p-v^  -  p-v)  =  0  (94) 


Let  r  represent  the  local  stagnation  temperature,  and  let  represent  the  local  static 

temperature.  Introducing  the  Perfect  Equation  of  State  and,  the  assumption  that  the  gas  is 
calorically  perfect.  Equation  (94)  becomes: 

- R(T-T)  -  —  +  R(T  -  T)  =  -^L--R(T-T)  -  —  =0  (^>5) 

(Y-1)  *  z  2  ^  ^  (Y-1)  ^  *  2 
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Under  the  assumed  gas  properties,  we  may  write: 


It  immediately  follows  from  Equation  (95)  that  at  every  point  in 


-l—‘\yRT  - 
(Y-1)  I  * 


Let  us  turn  now  to  the  Riemann  Invariant  on  characteristics  of  Type  I;  and  recall  that 
it  has  a  fixed  value  everywhere  in  S*,  determined  by  Equation  (93),  as  a  result  of  the  Type  II 
characteristics  being  straight  lines.  We  find  that: 


B(a„yC  + 
r  (Y“1)  ^ 


4  2 

-  C 

(Y-1)^ 


Hence  the  stagnation  temperature  T  may  be  written  in  terms  of  Riemann  Invariant 
BjiOjj)  and  local  sound  speed  C  at  every  point  in  Equations  (97)  and  (98)  imply  that: 


1  fY+ll  2 

T  =  —  •  (Y-D- -  -  2C-R/ap  H  ^  •C’ 

s  yR  2  ^  "  I  Y-l| 


Differentiating  this  expression  with  respect  to  sound  speed,  we  find  that: 


2  f  Y+1 1 

dc  yR  lY-lJ 


(100) 


Noting  Equations  (91),  (92),  and  (93),  we  find  that  —  =  0  at  plane  0-0  and 

dC 

decreases  monotonically  with  decreasing  C.  Hence,  Equation  (100)  implies  that  T* ,  the 
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supremum  of  T  ,  is  located  at  the  back  face  of  shock  ^  seen  in  Figure  7 ,  and  propagates 
with  it;  and  T  ^  ,  the  infimum  of  T  ,  is  permanently  located  at  the  plane  of  choked  flow. 
It  follows  from  Equations  (91),  (92),  (93),  and  (99)  that: 

T  =7’  (lOl.i) 

s*  0 


r* 


(101.ii) 


Equation  (99)  establishes  that  T  is  a  strict  function  of  C  on  Type  I  characteristics. 
Our  problem  is  to  calculate  the  mapping  of  T  on  domain  ®  in  the  (x,  t)  plane.  Consider 
Figure  6  again  and  take  account  that  in  the  case  of  the  backward-facing,  simple  wave,  the 
roles  of  the  and  F^  characteristics  are  reversed  from  those  in  the  mapping  portrayed 

there.  For  the  backward-facing  expansion  wave,  every  characteristic  F^  maps  to  the  same 

single  curve  R  in  the  (  p  ,U)  plane.  In  consequence,  ,  as  defined  in  Table  17,  is 

fixed  for  every  point  on  R .  Now  note  that  selection  of  a  point  p  in  ®  establishes  a  umque 

Oj  and  Ojj  on  a  characteristic  F^  ,  and  thus  establishes  a  point  p  on  R.  Clearly,  this 

is  a  unique  point  in  the  ( p  ,U)  plane;  hence  both  p  and  U  are  umquely  determined.  We  have 
simultaneously  uniquely  determined  C  ,  as  can  be  seen  by  utilization  of  the  fact  that 

B/o^)  is  a  constant  function  on  but  then  Equation  (99)  may  be  used  to  establish  a 

unique  value  of  T  which  corresponds  only  to  our  unique  value  of  C. 

Continuing,  the  unique  Oj  at  point  p  establishes  the  fixed  value  for 
corresponding  to  the  particular  F^  passing  through  p  ,  upon  which  U  and  C  are  fixed, 

since  a  Type  II  characteristic  is  a  straight  line  in  @.  Hence  /  1  F^  =  />  for  F^^  passing 
through  p  in  S'.  It  is  at  once  clear  that  selection  of  point  p  in  S  establishes  /  1  F^^  in 
its  entirety  for  F^  passing  through  point  p  . 

Letting  C  and  p  be,  respectively,  the  sound  speed  and  gas  density  at  the  plane  of 

choked  flow,  we  note  that  these  quantities  are  fixed  for  all  time.  But  the  isentropic  constraint 
implies  that  C  and  p  at  point  p  are  related  to  and  p^  according  to: 

2 

(102) 
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The  preceding  equation  shows  that  the  value  of  C  corresponding  to  point  p  on  JR  is  in 

one-one  correspondence  with  p  at  that  point.  Defining  T  as  the  mapping  of  domain  @  to 

the  set  of  (  T  ,  U)  pairs  *  corresponding  one-one  with  the  (  p  ,U)  pairs  comprising  curve  R , 
we  have  the  following  material  equivalence: 

/  r  =  T  r  im 


By  Equation  (102)  we  have  thus  shown  that  T  1  is  also  established  in  its  entirety  in 

0)  for  passing  through  point  p  .  Since  T  1“  F^  is  known  at  plane  0-0,  the  plane  of 

choked  flow,  and  is  stationary  there,  we  assert  that  the  properties  of  the  simple  -  wave 
mapping: 

f:  0  ^  R 


permit  the  exact  calculation  of  the  entire  stagnation  temperature  mapping  on  0. 

Note  that  this  wave  is  centered  at  point  (0,  0)  in  the  (x,  t)  plane,  and  let  us  refer  to 

this  point  as  0  .  This  means  that  every  characteristic  Fj^  in  domain  0  must  have  point 

0  as  a  limiting  point.  Inspection  of  Table  17  shows  that  AjjiOj)  =  0  everywhere  in  0, 

and  hence  is  known  a  priori  for  every  value  of  in  S*,  as  a  result. 

In  order  to  construct  the  mapping: 

T:0  ^  R 


we  need  to  know  the  stagnation  temperature,  ,  at  plane  0-0.  If  we  wish  to  construct  / 

in  conjunction  with  T  ,  the  stagnation  pressure,  ,  must  be  known  at  plane  0-0,  as  well. 

Given  these  quantities,  let  us  select  any  T  such  that  F  ^  ^  T  <  and  solve  the 
following  quadratic  equation  for  C: 


f  Y+1^  2 

J--  -c"  -  2C-B/o„)  *  (y-D- - yR  T  =  0 

Y“1  j  2  ® 


(104) 


Refer  to  Equation  (99). 
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The  root  for  which  —  ^  0  ,  as  determined  from  Equation  (100),  is  the  unique  value  for  C 
dC 

corresponding  to  T  for  every  (x,  t)  pair  in  &  mapping  to  T  .  Now  is 

determinate  from  Equation  (93);  hence  U  is  determinate  from  the  Riemann  Invariant  in  the 
first  row  of  Table  17.  This  establishes  the  appropriate  point  p  in  the  range  of  T  .  By 

Equation  (102),  point  p  in  the  range  of  /  is  also  established.  One  can  see  that  there 

exist  unique  U(p)  and  C(p)  such  that  there  is  a  family  2  of  lines  expressed 

as: 

=  |jC-(t/(^)-C(p))-f  = 


and  such  that  /  \  j  2  ~  ^  th^^e  exists  a  such  that  , 

where  F^^  is  the  unique  characteristic  of  Type  II  mapping  to  point  p  ,  for  some  member 
of  {1}  .In  our  case,  [Const]  =  AJo.)  =0  ;  and  since  F_  corresponding  to 

a  a=l,  2,  ...  a  u  I  u 

point  p  must  be  limited  at  6  ,  we  conclude  that  it  is  completely  described  in  domain  ®  by: 
X-[U^)-C(p)\t  =  0  (105) 


Thus,  when  T  has  been  established,  one  may  establish  a  point  p  in  the  range  of  function, 
T  ;  and  thence  establish  the  entire  mapping  T  t  Tjj  =  p  .  Since  we  are  free  to  pick  any 
r  such  that  r  ^  ^  ^  entire  mapping: 

T:Si  -*  R 


is  the  union  of  mappings  T  \  F^  ,  a  procedure  for  the  exact  calculation  of  the  mapping  of 
T  on  domain  S)  is  established. 

S 

It  is  not  necessary  for  the  calculation  of  the  exact  mapping,  T  ,  to  establish  the 
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trajectories  of  the 

characteristics.  However,  we 
may  wish  to  establish  these 
trajectories  in  the  (x,  t)  plane 
in  order  to  visualize  the 

properties  of  T  .  Consider 

point  o  in  Figure  8, 
and  let  us  suppose  that  J 
is  known  exactly  there.  This 
means  that  U  and  C  are 
exactly  known  there,  as  well, 

since  is  known  and 

fixed  at  every  point  in  S*  ^ 
Suppose  that  U  and  C  are 

known  for  another  F^ 
characteristic  in  domain  S*. 


X 

Figure  8 


Our  problem  is  to  calculate  the  coordinates  of  the  intersection  of  the  new  F^^  characteristic 


with  the  F^  characteristic  passing  through  point  o. 


Recalling  that  parametrizes  F^  ,  and  letting  A  o^  be  the  difference  between 


the  values  of  at  points  o  and  q,  we  note  that  the  Mean  Value  Theorem  guarantees  the 
existence  of  a  point  m  on  F^  between  o  and  q  such  that: 


x(q)  -  x(p)  =  [f7(w)  +  C(m)]-[f(^)  -  r(o)] 


(106) 


But  is  fixed  throughout  ®.  Hence,  Equation  (106)  becomes: 


x{q)  -  X(p)  = 


f  Y-3) 


•C(»i) 


\t{q)  -  f(o)] 


(107) 


Refer  to  Equations  (100)  and  (104). 
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We  shall  assiime  that  C  is  infinitely  differentiable  with  respect  to  .  In  this  case,  there 

exists  a  real  number,  K,  such  that  0  ^  K  ^  I  ,  and  such  that  the  following  Taylor’s 
expansions  hold: 


C'(m)  C"{jn)  2  2  C'\m)  3  ^  3 

C(o)  =  C(,m)  -  t  -  —^-X-’cAo,)’  * 


(108) 


and 


C(9)  =  C(m) 


C'(m) 


1! 


■(l-K)-Ao,  +  .^^•(l-X)"-(Aoy  * 


C^\m)  3  3 


(109) 


Letting 
Cip)  = 


C(o)  +  C(^) 


we  may  write: 


C(p)  =  C(m)  +  .^•(1-2J0  A<t,  .  .£|^-(l-2X+2x")-(Aa,)'  . 


C"'(m) 

2-3! 


(l-3«:+3A:^-2^^)*(Aa^)^  + 


(110) 


If  X(^)  is  the  predicted  value  of  x  at  point  q  using  C(p)  in  place  of  the  correct  mean 
value  for  C,  Equation  (107)  implies  that: 


Xiq)  =  x{o)  + 


U-ij 


•C(p) 


•[%)  -  t(p)\ 


(111) 
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Apparently,  Equations  (107)  and  (111)  imply  that: 


X{q)  -  x(q) 
C"'(m) 


^  Y-3^ 

i  Y-lj 

C\m)  ^  i  *  2 

(l-2iO'Aa,  +  -^'(l-2ii:-H2^:^-(Ao/ 


2-1! 


2*3! 


•(l-3A'+3ii:^-2A:V(Aa/  +  - 


•[%)  -«(o)] 


(112) 


Ao. 


We  note  that  ^t(q)  -  f(o)| 

(Aap{%)  -  t(p)]\  ^ 

Consequently,  the  following  inequality  holds 

|y-3 


.  This  inequality  implies  that: 


Ao 


xiq)  -  Xiq) 
d"{jn) 


2*3! 


Y-1 

(\-3K+3K^-lK\{LOjf 


C^(m)  2 

'•(1-2^0*(Ao/ 


2-1! 


C"{m) 


2-2! 


(i-2^:+2A’^- (A  op¬ 


ens) 


The  lead  term  of  Equation  (113)  shows  that  the  use  of  the  average  value  for  sound  speed  C 
on  an  interval  of  in  ®  results  in  a  second  order  error,  expressed  as  0(Aop  ,  in  Xiq). 

If  the  coordinates  of  point  o  in  @  are  known,  and  T  is  known  there;  and  T  can  be 
established  at  some  other  point  q  on  F^  passing  through  point  o,  then  the  coordinates  of  q 

in  can  be  established  to  0(A  .  Note  that  the  F^  characteristics  are  centered  at  the 

origin  of  coordinates.  Since  they  are  straight  Imes  in  we  may  write: 
xiq)  =  [f/(q)  -  Ciq)\tiq)  (114) 


Since  B^Ojj)  is  fixed  throughout  3),  Equation  (114),  which  is  exact,  becomes: 

(  Y+ll 


xiq)  = 


Y-lj 


'Ciq) 


tiq) 


(115) 
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Let  At  =  t(q)  -  tip)  and  Ax  =  xiq)  -  x(o)  .  Then,  to  within  0(Aop 
following  pair  of  equations  exists: 
xiq)  =  xio)  +  Ax 


tiq)  =  tio)  +  - 


Ax 


xiq) 


{  Y-3 

Y-1 


Cip) 


Y-1 


■Ciq) 


,  the 


(116.i) 

(116.ii) 


If  point  [xio),  tio))  is  established  on  a  Tj  characteristic,  then  Equations  (116.i)  and 


(116.ii)  establish  point  [xiq),  tiq))  on  that  same  characteristic  to  within  OiAOjf  . 


7.  RESULTS  FROM  EXAMPLE  CALCULATION  FOR  AN  INFINITE  RESERVOIR 
SUSTAINING  STAGNATION  PRESSURE  OF  65.0  PSIA  AND  STAGNATION 
TEMPERATURE  OF  740  “R  AT  PLANE  0-0 


Let  us  proceed  with  the  calculation  of  a  backward-facing  expansion  wave  issuing  into 
a  void  main  connected  to  an  infinite  reservoir  containing  standard  air,  and  maintaining  a 
stagnation  pressure  of  65.0  psia,  and  a  stagnation  temperature  of  740  ”  R  at  plane  0-0,  the 
plane  of  choked  flow.  A  characteristic  net  on  domain  S)  in  the  (x,  t)  plane  for  this  system, 
represented  by  Equations  (27.i)  and  (27.ii),  is  illustrated  in  Figure  9  on  the  following  page. 

The  initial  and  boundary  conditions  are  given  at  point  (0,  0)  in  3);  with  x  =  0  being  located 

at  the  plane  0-0.  The  line  (0,  t)  defines  the  lead  characteristic  in  The 

thermodynamic  state  is  constant  on  this  line  since  constant  stagnation  pressure  and 

temperature  are  perpetually  maintained  at  plane  0-0.  Hence,  /  =  0  throughout  3),  and  U 

and  C  are  each  fixed  on  each  F^  characteristic.  It  follows  that  each  F^^  characteristic 

is  a  straight  line  in  the  (x,  t)  plane;  but  each  F^  characteristic  carries  a  different 

thermodynamic  state,  and  has  a  different  slope.  The  wave  is  centered  at  (0,  0),  i.e.,  the 
characteristics  all  emanate  from  (0,  0).  Thus,  the  example  wave  is  a  centered,  backward- 
facing,  simple  rarefaction  wave.  The  F^  characteristics  span  the  wave,  and  fail  to  be 

straight  lines.  However,  they  may  be  calculated  to  within  0(Ao^)^  as  explained  in  Section  6. 
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Figure  9 

We  re-emphasize  that  it  is  not  necessary  to  calculate  the  trajectories  of  the 

characteristics  in  order  to  establish  function  /  ,  or  its  material  equivalent,  the  temperature 
mapping  T  . 

There  are  three  soUd  lines  illustrated  in  Figure  9  which  are  characteristics  of  type 
in  the  domain  of  the  example  wave.  Each  of  these  characteristics  is  marked  with  five 

points  from  domain  ®  of  the  wave.  The  points  associated  with  the  solid  circles  in  Figure  9 
are  designated  with  >  and  the  points  associated  with  the  squares  are  designated 

with  ,  •  •  .  ,  etc.,  where  the  index  increases  with  increasing  x  and  t.  When  the 
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shock  front  y  is  approached  in  domain  S*,  the  Tj  and  Tjj  characteristics  become 

coincident  in  the  limit.  The  solid  line  marked  with  the  inverted  triangles  represents  that  limit. 
Now  the  range  of  /  on  for  the  example  wave,  is  the  single  curve,  R  ,  in  the 
(p,  U)  plane  -  illustrated  in  Figure  10.  Considering  the  developments  of  Section  5,  it  is 

clear  that  /  has  no  inverse.  All  of  the  points  on  a  given  characteristic  map  to  a  single  point 

on  the  curve.  For  example,  all  five  domain  points  indicated  by  solid  circles  in  Figure  9  map 

to  the  point  associated  with  the  single  solid  circle,  designated  as  n  ,  on  the  curve  in  Figure 

10.  The  lead  characteristic  coincides  with  the  ordinate  passing  through  (0,  0),  and  is  not 
marked  in  Figure  9.  However,  the  range  point  to  which  it  maps  is  illustrated  as  the  hollow 
circle  labelled  0  in  Figure  10. 


Figure  10 
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We  observe  that  escape  velocity,  17^  ,  discussed  on  page  42,  is  proportional  to 
at  plane  0-0,  given  by  Equation  (91).  The  relationship  may  be  written  as: 

Y+l'^ 


=  Bfo^)  = 


Y-lj 


‘U 


cr 


(117) 


In  Figure  10,  the  inverted  triangle  denotes  the  limit  of  characteristics.  One  can  deduce  from 
its  position  that  in  our  example  wave,  =  7303.4726  ftjs&c  while  =  1217.33  ^/sec 

The  stated  magnitudes  of  these  velocities  are  in  precisely  the  ratio  given  by  Equation  (117) 
when  Y  for  standard  air  is  incorporated. 

It  is  clear  that  the  limit  velocity  at  plane  0-0  is  substantially  less  that  escape  velocity 
for  3).  The  equation  for  limit  velocity  at  plane  0-0  in  terms  of  both  BjiPjj)  and  is 


given  as  follows: 


U  = 

10  M 


W-1 


{  Y+ll 


Y-1 


•u 


(118) 


Incorporatmg  y  for  standard  air,  we  find  that  =  2981.84  ftjsitc  in  our  example  wave. 

One  can  see  that  is  larger  than  by  a  factor  of  2.4495  -  the  effect  being  due  to  the 

flow  work  done  to  form  the  shock  at  t/at  the  time  of  initiation  of  the  expansion  wave. 

Noting  that  /  symbolizes  the  subject  wave  function,  we  have  shown  that  imder  the 
hypotheses  laid  down  in  Section  3,  the  representation  of  the  mapping: 

f  :2}  ^  R 


given  from  characteristic  theory  is  exact.  We  have  also  shown  that  on  each  Tjj 
characteristic,  /  1“  =  T  I  .  This  means  that  if  one  were  to  select  a  point  in  such 

as,  {/>,}, -,1  from  Figure  9,  the  unique  pair,  (T  ,  17) ,  which  corresponds  to  it  may  be 
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calculated  precisely  and  directly.  Since  the  characteristics  are  straight  and  centered  at 

(0,  0),  there  exists  a  constant,  A  ,  such  that  0  ^  ^4  <  BjiOjj)  ,  and  such  that  the 
following  equation  pair  holds  in 

-  =  -C  ^  U=A 
t 

—■C  *  U  =  Bio„)  (»»•«) 

Y-1  ^  ^ 


Now  B^o„)  is  constant  throughout  and  determinate  from  the  given  stagnation  tem¬ 
perature  at  plane  0-0,  using  Equation  (93)  from  Section  6.  For  the  selection  of  a  point  in  Q), 

X 

such  as  ,  it  is  clear  that  A  is  determined  by  -  for  that  point.  Hence  C  and 

U  are  each  uniquely  determined  by  the  following  equation  pair: 


(120.i) 

(120.ii) 


The  stagnation  temperature,  T  ,  is  immediately  determinate  by  means  of  Equation  (99) 
from  Section  6,  as  a  result. 

Under  the  mapping,  T  ,  there  is  a  unique  pair,  (T  ,  17) ,  corresponding  to  each 
choice  of  (x,  t)  in  Hence  T  is  a  function  of  x  and  t,  and  our  wave  may  be  represented 

in  this  fashion.  Note  that  in  view  of  the  material  equivalence,  /  [■  F^  =  T  F^  ,  and  the 
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condition  that  J  =  0  ,  every  point  on  the  characteristic  passing  through  a  given  point 
in  %  maps  to  the  same  value  of  T  .  If  our  choice  of  (x,  t)  had  been  point  in 

Figure  9,  for  instance,  then  the  points,  tP,},--i^34^  >  represented  by  the  solid  squares,  as 
well  as  any  of  the  other  points  on  the  F^^  characteristic  passing  through  »  would 

map  to  the  same  value,  T  =  2820®i?  . 

The  wave  function  solving  our  example  problem  is  represented  as  T  (x,  t)  in  Figure 

11,  seen  on  the  following  page.  Each  contour  label  is  the  distance,  x,  in  feet,  from  the  plane 
of  choked  flow,  i.e.,  plane  0-0.  Again,  we  note  that  this  representation  is  exact.  Hence, 

Figmre  11  is  the  representation  of  the  function  T(t)  at  each  of  the  indicated  distances 

from  plane  0-0.  Stagnation  temperature,  T  ,  depending  upon  in  Figure  9,  is 

the  set  of  points  represented  by  the  solid  squares  in  Figure  11,  which  are  sequenced  from  left 

to  right  with  increasing  index  i.  The  common  value,  T  =  2820®/?  ,  which  applies  at  each  of 

the  points,  lP,},=i^34^  >  appears  on  each  of  the  five  contours  to  which  each  of  these  points 
corresponds. 

Referring  to  Equation  (119.i),  one  can  see  that  constant  A  is  the  propagation  velocity 
of  a  particular  point  in  the  range  of  the  wave  carrying  a  particular  value  of  stagnation 
temperature  T  .  In  our  example  case,  we  directly  calculate  that  A  =  5476.8592  ftlsec  for 

the  Fjy  characteristic  marked  by  points  {F,},=i^34^  in  Figure  9.  Employing  Equations 

(120.i)  and  (120.ii),  and  noting  that  -  7303.4726  ^/sec  ,  we  find  that 

C  =  304.43557  ftjs&c  and  U  =  5781.29468  ^/sec  .  Turning  to  Equation  (99),  we 

calculate  that  T  =  2820.684®/?  ,  which  is  sufficiently  close  to  the  previously  determined 
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value  of  2820°/?  .  This  calculation  of  course,  can  be  performed  at  any  point  in  Q). 


It  is  clear  that  each  value  of  T  in  the  range  of  our  wave  has  a  different  velocity  of 
propagation;  the  range  of  velocities  being  0  ^  .  The  value  of  T  at  the  critical 

flow  point,  i.e.,  plane  0-0,  is  given  by  Equation  (lOl.i).  At  this  location,  T  =  T,  . 
However,  the  value  of  T  at  the  rear  face  of  the  shock  is  6  *  in  standard  air 
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according  to  Equation  (lOl.ii).  Here,  we  write:  =  T*  .  The  extreme  temperatures  in  the 

range  of  our  wave,  as  well  as  that  given  for  our  example  point  p  ,  along  with  their 
velocities  of  propagation  are  given  in  the  following  table: 


T 

S 

Velocity  of  Propagation 

Shock 

4440 

7303.47260  fl/sec 

Point  p 

2820  "R 

5476.8592  ft/sec 

Critical  Point 

740  °R 

0.0  fVsec 

It  is  at  once  apparent  that  the  critical  point  in  the  range  of  the  wave,  represented  by  the  last 
line  of  the  table,  has  no  propagation  velocity  at  all.  It  does  not  advance  into  the  infinite  main. 
It  is  on  the  lead  characteristic  of  the  wave,  and  is  permanently  located  at  plane  0-0  from 
initial  time  forward. 

One  can  see  by  examination  of  Figure  11  that  the  shape  of  T  (f)  is  different  for 

each  of  the  five  linear  distances  from  plane  0-0  indicated  by  the  contour  labels.  This  is  the 
result  of  each  part  of  the  wave  having  a  different  propagation  velocity.  The  distribution  of 

r  (r)  ,  at  each  distance,  is  an  artifact  of  the  initial  flow  work  done  to  create  the  shock 

which  proceeds  into  the  vacuum  of  the  main.  The  vertical  lines  mark  the  time  of  appearance 
of  the  shock,  following  initiation  of  the  wave  at  time  zero,  at  each  of  the  five  indicated 
distances  from  plane  0-0. 

Let  us  turn,  now,  to  the  static  pressure  function,  p{x,t)  ■  At  each  distance  x 

from  plane  0-0,  /»  =  0  at  initial  time,  and  asymptotically  approaches  a  limit  given  by  the 

following  formula: 


Pcr^Po 


2  ^ 


u 


cr 


2 


Y 

Y-1 


(121) 


where  and  are  precisely  given  by  Equations  (91)  and  (118),  respectively.  This 
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limit  exists  due  to  being  the  lower  limit  of  velocity  for  domain  S).  For  our  example 

case,  we  calculate  =  34.3381  psia  from  Equation  (121). 

Figure  12  shows  j?  (f )  at  each  of  the  same  five  distances  from  plane  0-0  given  by 

the  contour  labels  in  Figure  11.  Those  distances  are  represented  by  the  contour  labels  in 
Figure  12,  just  as  they  were  in  Figure  11.  The  indicated  points  on  these  contours  are  the 
images  of  the  similarly  marked  points  on  characteristics  shown  in  Figure  9;  as  was  the  case 

for  the  r  points  in  the  range  of  T(t)  shown  in  Figure  11.  The  grad  symbol  was  made 

hollow  in  this  plot  so  that  the  image  points,  3  >  could  be  represented  despite  the 

image  positions  being  so  compacted  due  to  the  relatively  gradual  rate  of  rise  of  static 


Figure  12 
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pressure.  The  static  pressure  limit,  i.e.,  critical  static  pressure  is  denoted  by  <E>  . 

On  the  interval  over  which  the  analysis  was  performed,  Figure  12  clearly  shows  the 

asymptotic  behavior  of  /?(?)  at  one  foot  distance  from  plane  0-0;  although,  even  at  the  one- 

foot  distance,  it  has  only  achieved  91.7%  of  its  limit  value  after  10  milliseconds. 

The  variation  in  velocity  of  propagation  throughout  @  has  considerable  effect  on  the 

form  of  (f )  as  distance  from  plane  0-0  increases.  The  form  of  (r)  at  25  feet  from 

plane  0-0  is  dissimilar  to  its  form  at  1  foot  distance  -  its  approach  to  the  limit  being  much 
more  gradual  at  25  feet,  as  can  be  seen  from  the  lower  right-hand  contour  in  Figure  12. 

Figure  13  represents  T  (x,f)  as  T  (jr)  at  various  times  after  initiation  of  the 
wave  at  plane  0-0.  In  this  figure,  the  contour  labels  represent  the  time  in  seconds  after  the 


Figure  13 
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initiation  of  the  expansion  wave,  while  the  contours,  themselves,  represent  the  distnbution  of 
stagnation  temperature.  In  this  figure,  the  vertical  lines  give  the  shock  position  at  each  of  the 
times  indicated  by  the  contour  labels.  The  solid  squares  indicate  the  stagnation  temperature 
immediately  behind  the  shock  at  each  of  these  times. 

The  stagnation  temperatme  immediately  behind  the  shock  retains  the  identical  value  of 
4440°  R  at  each  of  the  shock  positions  shown  in  Figure  13.  However,  it  is  clear  that  the 
stagnation  temperature  decreases  at  each  location  after  the  passage  of  the  shock,  eventually 
approaching  the  lower  limit  of  740°  R  marked  by  the  dotted  line.  At  each  distance,  x,  and  at 
this  limit,  the  flow  velocity  is  precisely  critical. 

Finally,  the  stagnation  temperature  in  the  critical  flow  existing  at  plane  0-0  is  the 
result  of  the  work  of  expansion  to  accelerate  the  fluid  from  its  quiescent  condition  in  the 
reservoir  to  its  state  at  plane  0-0.  This  acceleration  is  accomplished  by  an  isentropic 
expansion  wave  which  propagates  into  the  subsonic  flow  in  the  reservoir,  in  the  opposite 
direction  from  the  rarefaction  wave  in  the  main.  The  result  is  that  the  stagnation  temperature 
and  pressure  at  plane  0-0  are  each  always  less  than  the  stagnation  temperature  and  pressure  in 
the  reservoir.  For  our  reservoir,  which  extends  from  the  origin  to  ~  on  the  x-axis,  and  is 
infinite  in  extent  in  every  direction  perpendicular  to  the  x-axis,  and  in  which  the  flow  is  one¬ 
dimensional,  this  expansion  wave  is  confined  to  plane  0-0  and,  in  this  limiting  case,  becomes 
a  stationary,  isentropic  expansion  shock  *.  One  can  deduce  that  the  Second  Law  of 
Thermodynamics  is  not  contradicted  in  this  instance 

The  stagnation  pressure  and  temperature  at  plane  0-0  are  maintained,  respectively,  at 

=  65.0  psia  and  =  740.0®R  after  initial  time  in  our  example  problem.  Under  this 
condition,  the  corresponding  static  pressure  and  static  temperature  at  plane  0-0  are  each 
=  34.338  psia  and  T  =  616.667°!?  .  The  expansion  wave,  which  must  be  present  at 
plane  0-0,  requires  that  the  stagnation  pressure  and  temperature  in  the  reservoir  remain  at 
p  =  82.419  psia  and  T  =  791.940°!?  ,  respectively,  for  f  ^  0  ,  in  order  to 

maintain  the  specified  values  of  p^  and  at  plane  0-0. 

8.  CONCLUSIONS 


A  one-dimensional  analysis  of  the  flow  resulting  from  the  expansion  of  a  calorically 
perfect  gas  from  an  infinite  reservoir  into  a  void  main  implies  that  the  rarefaction  wave 
associated  with  this  expansion  is  preceded  by  a  shock  propagating  into  the  vacuum  of  the 
main  at  a  rate  given  by  Equation  (1 17)  on  page  55.  The  argument  leading  to  this  conclusion 
is  found  in  Section  5  on  pages  40-42.  This  argument  also  implies  that  the  velocity  of 
propagation  of  the  shock  is  the  least  upper  bound  of  the  velocity  of  propagation  of  a 
disturbance  in  the  domain  of  the  rarefaction  wave.  If  the  specific  heat  ratio  for  standard  air  is 
applied,  then  the  propagation  velocity  of  the  shock  is  six  times  the  critical  velocity 
corresponding  to  the  thermodynamic  state  at  plane  0-0,  seen  in  Figure  1,  Section  3. 

‘  The  development  of  this  model  is  given  in  the  Appendix. 

^  See  the  development  of  Equations  (Al),  (A2),  and  (A3)  in 
the  Appendix,  as  a  limiting  case.  62 


Concomitantly,  the  stagnation  temperature  immediately  following  the  shock  is  given  by 
Equation  (101. ii).  The  argument  leading  to  this  equation  is  given  in  Section  6,  pages  42-46.  If 
the  specific  heat  ratio  for  standard  air  is  applied,  the  stagnation  temperature  following  the 
shock  is  found  to  be  six  times  the  stagnation  temperature  at  plane  0-0,  as  well. 

The  roles  of  the  First  and  Second  Laws  of  Thermodynamics  in  the  development  of  the 
transient  flow  characterized  by  this  wave  are  developed  in  Section  3,  and  lead  to  Equation 
(95)  on  page  44;  and  to  Equations  (99)  and  (100)  on  page  45  in  Section  6.  These  equations, 
in  turn,  lead  to  the  conclusion  that  the  local  stagnation  temperature  in  the  domain  of  the 
rarefaction  wave  is  strictly  dependent  upon  local  sound  speed;  and  increases  as  sound  speed 
decreases,  reaching  its  limit  at  the  back  face  of  the  shock.  This  limit  temperature  does  not 
change  over  time,  and  appears  immediately  following  the  shock  at  each  downstream  location 
in  the  main. 

Apparently,  the  heat  of  compression,  i.e.,  the  flow  work  done  to  form  the  shock  when 
the  flow  is  initiated  at  plane  0-0,  produces  the  thermal  spike  which  one  calculates  fi’om  the 
application  of  characteristic  theory  to  the  flow  interior  to  the  rarefaction  wave.  Since  the  flow 
is  inviscid  in  the  analysis  presented  in  this  report,  the  thermodynamic  process  is  reversible 
everywhere  in  the  interior  of  the  wave;  hence,  the  calculated  spike  is  not  the  result  of 
frictional  heating  in  the  established  flow.  Even  though  the  one  dimensional  and  inviscid 
hypotheses  represent  substantial  simplification  of  the  actual  physical  flow  found  in  the 
transient  blow-down  facility,  they  embody  sufficient  physics  to  imply  that  the  thermal  spike 
observed  there  is  an  artifact  of  the  flow  work  done  to  establish  the  shock. 

The  one-dimensional  hypothesis  implies  that  the  flow  in  the  reservoir  must  be 
quiescent,  i.e.  in  a  state  of  rest,  and  therefore  steady,  even  at  the  initiation  of  the  rarefaction 
wave.  This,  in  fact,  is  not  true,  even  in  an  infinite  reservoir.  The  flow  in  the  neighborhood  of 
the  inlet  to  the  main  must  have  some  radial  component,  and  must  be  established  by  an 
expansion  wave  which  proceeds  into  the  reservoir,  however  diffuse  that  wave  must  be.  This 
means  that  the  stagnation  temperature  at  plane  0-0,  where  the  flow  is  critical,  is  actually  less 
than  the  stagnation  temperature  in  the  quiescent  state  in  the  reservoir. 

By  insisting  that  the  flow  be  one-dimensional  in  the  infinite  reservoir,  one  must  admit 
the  existence  of  an  isentropic  expansion  shock  at  plane  0-0.  It  is  shown  in  the  appendix  that 
the  isentropic  shock  does  not  contradict  the  Second  Law  of  Thermodynamics  under  the 
hypotheses  of  this  problem;  rather,  it  exists  as  the  limiting  case  of  the  expansion  wave  which 
proceeds  into  the  reservoir.  Hence,  the  representation  of  Ae  physics  of  the  development  of  the 
backward-facing  expansion  wave  is  not  fundamentally  altered  by  the  one-dimensional 
hypothesis.  The  application  of  characteristic  theory  to  the  two-dimensional  equations  of 
motion  would  allow  one  to  characterize  the  expansion  wave  which  proceeds  into  the  reservoir, 
and  therefore  estimate  its  true  effect. 
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APPENDIX 


Consider  again  the  adiabatic,  infinite  reservoir  and  void  main  of  Figure  1.  The 

reservoir  has  infinite  extent  in  .  Reservoir 

the  horizontal  direction  to  the  I/*  Quiescent  State 

left  of  plane  0-0,  and  in  any 

radial  direction  relative  to  its  / 1 

axis.  Let  us  extend  a  4  \ 

cylindrical  control  volume,  <2’,  |  pPiane  of  Choked  Flow 

of  the  same  diameter  as  the  |/  Shock 

main  from  plane  0-0  to  any  0  . 

distance,  L,  to  the  left  of  . — . 

’  ’  ( I - = - >  •  Vac 

plane  0-0.  We  suppose  that 

the  ri^t  face  of  ^  is  pw._sv.— 

coincident  with  plane  0-0,  |  V. - - - - ^ 

and  that  ^  is  coaxial  with  the  |  Rarefaction 

main.  Let  the  following  ^ 

notation  apply:  a  rep-  | 


Shock 

Front 


Vacuum 


Rarefaction 


resents  the  surface  of  £ 
coinciding  with  plane  0-0, 
represents  the 

remainder  of  the  surface  of 


Gas  Expansion  From  Infinite  Reservoir  Into 
Vacuous  Main 


when 


is  excluded,  and  O  represents  the  entire  surface  of  2.  Let 


represent  the 


length  of  as  L  is  extended  to  -00.  Note  that  L  ^  does  not  represent  a  finite  number,  nor  an 

attainable  length.  Finally,  let  subscript  RO  apply  to  the  thermodynamic  state  in  the  infinite 
reservoir,  and  let  subscript  1  apply  to  the  thermodynamic  state  at  plane  0-0. 

Now,  the  thermodynamic  state  is  fixed  in  time  throughout  <2',  and  accumulated  mass  of 

gas  is  infinite  when  h=L^ .  Subsequently,  if  the  flow  is  one-dimensional,  the  following 
conservation  relations  hold  for  the  half-open  time  interval  0  <  t  <  °° 

Momentum 


f  f  [V(pV’^)  +  P'n]da  =  p 


p  +  p  U 

'^1  cr 


•O  =  0 
0-0 
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Mass  Flow 


J|  pV’fida  =  JJ^  pV'nda*  +  =  0 

L‘L__-  |k(0.)  =  0  A  p(o.)  =  p^}  for  P  >  0 

Entropy 


(A2) 


fUpnida  =:  s^-fl^^pv-tido. 


s-pU 


a 

cr  0-0 


=  0  (A3) 


L  =  L 


s  =  s 
1  so 


Equations  (A2)  and  (A3)  imply  that  when  L=L  ^  ,  the  steady  effluxes  of  mass  and 

entropy  from  are  balanced  by  the  integrals  of  vanishingly  small  fluxes  of  these  quantities 
over  surface  a  ^  of  S.  The  entrained  mass  and  entropy  in  are  each  infinite,  and  remain 

steady  over  the  specified  time  interval.  The  momentum  entrained  in  i?’  is  zero  over  the 
specified  time  interval  in  this  case.  Equation  (A3)  shows  that  the  flow  entering  the  head  of 

the  rarefaction  wave  at  plane  0-0  is  carrying  entropy  for  o  ^  f  <  «»  ;  and  hence  that 

s  =  throughout  domains  SS  and  and  their  intersection,  which  is  plane  0-0. 

Let  us  apply  the  First  Law  of  Thermodynamics  to  the  expansion  from  the  infinite 
reservoir  into  the  main.  Recall  that  h  defines  total,  or  stagnation  enthalpy.  The  flow  is 
quiescent  in  the  reservoir,  and  time-steady  at  plane  0-0  for  0  •  Now  the  total 

enthalpy  efflux  through  the  surface,  a ,  of  written  as:  J  j  (pV'Ayhda  ,  is  the  sum  of  two 

components;  these  components  being  expressed  as  Equations  (A4.i)  and  (A4.ii): 

(pV'Ayhda  =  JJ  (pVAyeda  +  jj  (pV'Ay(pv)da 
//  (py-^i‘do  - 

O 


(A4.i) 

(A4.ii) 
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Suppose  that  L  approaches  L  .  The  infinite  reservoir  is  in  a  quiescent  state,  so; 
f  f  f  De  CP  - 

-aT  ■  ■  lfjpv-«)-ed„  -  « 

L  =  L  =>  ivia  )  =  0  A  p(o  )  =  p„„  A  c(a  )  =  e  \for  peU  ‘O  >  0 

-«  I  ^RO  ^  R0\  ^  \  \  cr  0-0 

Consider  a  point  p  in  ,2’.  In  view  of  the  quiescent  state  of  the  reservoir: 

De 

pes^—(p)=0 


It  follows  that  when  L=L  ,  we  should  write: 

—OP  ^ 


It  is  apparent  that  when  L=L  ^  ,  the  efflux  of  total  energy  at  plane  0-0  no  longer  materially 

changes  the  total  energy  of  the  gas  contained  in  the  reservoir.  Additionally,  Equation  (A2) 
implies: 

ff  pV-nda  =  -p  f/  (A8) 

JJo+  1  cr  0-0 


Consequently,  Equations  (A4.i)  an  (A4.ii)  may  be  written  in  the  following  form: 
ffjpmhda  =  - 

If  (.P^mdo  =  pV^e^-o,^  *  P,‘'„’Pi‘''rVo 


(A9.i) 

(A9.ii) 


Equations  (A9.i)  and  (A9.ii)  imply  that  when  L=L  ^  ,  the  integral  of  the  stagnation  enthalpy 


over  the  surface  of  <2’ becomes: 

J/ (pP-^-fcdo  .  P,l'„-(p,-v,-p„-v^)-o„^  <  0 


(AID) 
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The  right-hand  side  of  Inequality  (A  10)  is  clearly  recognizable  as  the  flow  work  rate  which 
must  be  supplied  by  the  reservoir  to  accelerate  the  flow  from  its  quiescent  state  interior  to  the 
reservoir  to  its  critical  state  at  plane  0-0.  This  term  must  invariably  be  strictly  less  than  zero. 
Evidently,  the  First  Law  Thermodynamics  does  not  require  the  total  enthalpy  efflux  from  2  to 
vanish  in  case  the  reservoir  is  infinite  in  extent,  even  when  the  flow  at  plane  0-0  is  fixed  in 
time. 

It  is  clear  that  a  jump  discontinuity  in  the  total  enthalpy  at  plane  0-0  is  permitted, 
since  the  infinite  reservoir  maintains  an  inert  state  of  total  energy  in  spite  of  the  facts  that  it 
is  adiabatic,  the  flow  conditions  are  time-steady,  there  is  no  influx  of  total  energy  into  the 
reservoir;  but  there  is  an  efflux  of  total  energy  at  plane  0-0.  An  expression  admitting  a  jump 
of  finite  magnitude  is  given  by  Inequality  (A  10).  However,  the  magnitude  of  the  jump  is  not 
determined  by  this  inequality.  Rather,  given  that  the  thermodynamic  process  is  reversible,  it  is 
determined  by  the  entropy  conservation  relation  of  Equation(A3).  One  can  see  that  the 
infinite,  one-dimensional  reservoir  permits  an  isentropic  expansion  shock  at  plane  0-0  which 
satisfies  both  the  First  and  Second  Laws  of  Thermodynamics. 

Next,  we  consider  the  isentropic  process  at  plane  0-0.  Let  the  following  definitions 

hold: 


r  = 


(All) 


(A12) 


In  view  of  the  application  of  the  Perfect  Equation  of  State,  we  may  write: 

_ 

Pi  jci 
1 


Hence,  Equation  (Al)  implies  that: 


KO 


/ 


Pi 


\ 


=  0 


(A13) 


(A14) 


and  we  may  write: 


r 


1  + 


CT 


RT 

1 


(A15) 
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Continuing,  let  US  compare  with  the  critical  temperature  ratio  at  plane  0-0.  it  is 
apparent  that  the  isentropic  condition  on  domains  SE  and  as  well  as  their  intersection, 
implies  that  one  may  define  a  function  P(z)  such  that: 


P(z)  =  (z) 


Y-l 

Y 


(A16) 


where  z  is  an  arbitrary  pressure  ratio  which  is  infinitely  differentiable.  Now,  by  Taylor’s 
Theorem,  there  must  be  a  real  nxunber  (  between  1  and  r  such  that  the  following  relation 
must  hold: 


T  1!  2!  3! 


When  Equation  (A15)  is  introduced.  Equation  (A17)  becomes: 


^  ^  (Y-l) 

T 
1 


(Y-1).(Y+1) 


(Y-l) 


/  ,  N2 

u 


1!  yRT^ 


3! 


2y+1 

(  TT^ 

. . 

u 

Y  . 

cr 

yRT 

V 

2! 

3 


Yi?r 


1  <  C  <  ^ 


(A17) 


(A18) 


For  critical  flow  at  plane  0-0,  it  must  be  that: 

r2 


=  1 


Yi?r 

In  this  case.  Equation  (A18)  reduces  to: 


_  (Y+1)  ^  (Y-1)-(Y+1) 

r  "  2  "  6 


/  1  \ 


cj 


2y+1 


1  <  C  <  (Y+l) 


(A19) 


(A20) 
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Equation  (A20)  shows  that  ratio  of  reservoir  temperature  to  the  temperature  at  plane  0-0 
differs  from  critical  temperature  ratio  only  by  the  remainder  of  the  Taylor’s  Expansion  after 
the  first  three  terms.  Re-introducing  the  notation  of  Equation  (A12),  we  may  write  an 
inequality  which  expresses  the  difference  between  the  reservoir  temperature  ratio  and  critical 


(Y-IHY+1)  (A21) 

6 


temperature  ratio  at  plane  0-0.  This  inequality  is: 

2y+1 


(Y-1)'(Y+1) 


1 


Y+IJ 


<  P.  - 


(Y+1) 


Clearly,  is  steady  for  0  <  f  <  ~  ,  and  C  is  a  unique  number.  It  follows  from 


Equation  (A20)  that  is  fixed  for  0  ^  f  <  «»  as  well. 

Note  that  if  symbolizes  the  stagnation  temperature  at  plane  0-0,  and 

symbolizes  the  static  temperature  there,  when  the  flow  velocity  is  critical,  then: 


Y+1 


(A22) 


Equation  (A22)  is  the  aforementioned  critical  temperature  ratio.  Obviously,  is  fixed,  and 

T  >  T  for  0  ^  f  <  oo  •  If  the  working  fluid  is  standard  air,  then  the  critical 
no  0 

temperature  ratio  at  plane  0-0  may  be  calculated  as: 

=  1.200  (AIJ) 

T  2 
1 

However,  Inequality  (A21)  implies  that: 

0.01486  < 


Y-^1) 


<  0.16000 


(A24) 


Clearly,  the  difference  is  always  greater  than  zero. 

The  values  of  reservoir  stagnation  pressure  and  temperature  which  must  be  maintained 
upstream  of  the  expansion  shock  in  order  to  sustain  critical  flow  at  plane  0-0  can  be  directly 
calculated  from  the  conditions  imposed  by  Equations  (Al),  (A2),  and  (A3).  They  are  given  by 
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Equations  (A25)  and  (A26): 


(A25) 


(A26) 


